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Summary 


An approximate analysis is made of the sound field in a rectangular box excited by a spherical 
source. Two different kinds of characteristic frequencies are observed, one of which is the set of 
eigentones. A general picture is obtained of the relative importance of the various modes. High 
frequency eigentones are relatively unimportant. 


Sommaire 


On fait ’é6tude approchée du champ sonore excité dans une enceinte rectangulaire par une 
source sphérique. On a observé deux types différents de fréquences caractéristiques, dont l’un 
est constitué par la série des sons propres. On a établi un schéma général de importance relative 
des différents modes. Les sons propres de fréquence élevée sont relativement peu importants. 


Zusammenfassung 


Das Schallfeld, das ein Kugelstrahler in einem rechteckigen Raum erzeugt, wird naherungs- 
weise analysiert. Zwei verschiedene Arten von charakteristischen Frequenzen werden beobachtet. 
Eine von ihnen ist das System der Eigenténe. Ein allgemeiner Uberblick iiber die relative Bedeu- 
tung der verschiedenen Schwingungsformen wird erhalten. Hohe Eigenténe sind ziemlich unbe- 


No. 


deutend. 


1. Introduction 


If a continuous sound is made in the open air, 
sound waves travel radially outwards in every di- 
rection and the sound pressure at any point is in- 
versely proportional to the distance from the source. 
In the absence of absorption or reflection the sound 
field ultimately fills all space. 

If, however, there is a plane reflecting surface near 
the source all the sound which would have reach- 
ed the far side of the reflecting plane appears on the 
near side superimposed on the original sound field: 
The result is equivalent to two sound sources — the 
original and its image in the reflecting surface — 
and the sound pressures are not arithmetically 
additive but combine in accordance with their rela- 
tive phases to give a system of nodes and anti- 
nodes dependent on the frequency. 

Therefore the sound field has been folded back, 
so to speak, into half the space giving twice the 
average pressure but with a node and antinode 
pattern due to interference. 

If now the source is completely enclosed by re- 
flecting walls, the sound which would have filled 
a large space is reflected into the enclosed volume, 


and the condition is equivalent to an infinite num- 
ber of sources consisting of the original source and 
its images in the walls and the images of the images 
and so on. 

In any studio, the increase in loudness over open 
air conditions is due to this reinforcement and the 
loudness varies considerably with frequency and 
position. The acoustical properties of a studio do 
not depend on the mere presence or absence of 
standing wave patterns but on the nature of the 
patterns. A studio with small standing wave field 
is a dead studio. 

In this paper, a study is made of the sound re- 
ponse of a rectangular room in the first instance to 
a point source of impulsive sound. The advantages 
of the impulsive method are discussed by vAN DER 
Pot [1] and here it is sufficient to say that the 
spectrum of the impulse is similar to the spectrum 
of random or “white” noise. In “white”? noise no 
specific pitch predominates and the spectrum of 
the impulse is the same except for the relative phase 
of the component frequencies. 

In studio acoustics, the use of an impulsive sound 
is useful for rapid diagnosis. An instrument in which 
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a small paper diaphragm is ruptured by air pressure 
by means of an air pump is a convenient source. If 
such an “‘acoustic impulse generator” is operated in 
a completely dead room, the resultant sound heard 
is almost zero. In reverberant rooms, the response 
heard is characterised by the surroundings and re- 
sembles a pistol shot. 

Such tests demonstrate the importance of the 
auditorium as an essential part of programme mak- 
ing apparatus. 


2. General discussion 


2.0. General method 


In this section the sound pattern in a rectangular 
room is discussed in general terms. Section 3 is in 
essence a restatement of this section in mathematical 
terminology. 

2.1.0. General 

From a loudspeaker any desired signal can be 
sent into a room and the resultant sound field 
observed by a microphone. To simplify the issue 
the loudspeaker is taken as a small omnidirectional 
sound pressure generator and the microphone as an 
omnidirectional pressure microphone. 

In practice, a loudspeaker is much more nearly a 
velocity source, or in other words, the cone move- 
ment is not much affected by sound pressure. 
Pressure and velocity are proportional in a plane 
sound wave but not in the divergent wave from a 
small source and it is for this reason that a velocity 
microphone should not be used at too short ranges. 
It is, however, easy to calculate the velocity input 
to give any required pressure wave [2]. 

The loudspeaker-room-microphone combination 
may be treated as a three terminalelectrical net- 
work. By placing the microphone and loudspeaker 
in various positions in the room the whole res- 
ponse of the room can be determined. There are 
several methods of observing an unknown network 
and each method gives all the information required, 
that is, completely determines the network [1]. The 
best known method is the measurement of the 
steady state characteristic at allrelevant frequencies 
but the method preferred here is the impulse re- 
sponse method. 


2.1.1. Pulse response and echo time dia- 
gram 

In this method the magnitude and epoch of all 
the echoes of an impulse are observed experiment- 
ally by means of an oscillograph with a linear time 
base. The exact shape of the pulse is not important 
experimentally provided it is short enough to enable 
the separate echoes to be distinguished. From this 


echo time diagram can be deduced any response, ~ 
whether transient or steady state, as will be dis- — 
cussed in the next paragraph. 

An echo time diagram is shown schematically in~ 
Fig. 1. In this diagram A is the start of the time base _ 
and.at this instant the impulse of sound is trans- 
mitted. At B is the trace of the sound travelling © 
direct to the microphone. At C is the first echo and 
D, E, etc. are later echoes. The length AB is pro- 
portional to the distance between the loudspeaker 
and the microphone and the distance AC is pro- 
portional to the distance between the nearest image 
and the microphone and so on. 


P 


pulse pressure 
S 


nh 


Bama is 
time 
Fig. 1. Echo time diagram. 


The ordinates of the curve PQR are inversely 
proportional to the time after the impulse is sent. 
Strictly speaking, the time is measured from the 
time of departure of the impulse from the equi- 
valent point source, not the time of departure from 
the actual speaker diaphragm. The magnitude of 
each echo decreases inversely with time or distance 
as well as by air attenuation and reflection loss 
so that the echo pulses all lie below the curve PQR. 

Suppose now that a continuous tone is sent from 
the speaker at constant level and starting abruptly 
at zero time. The resulting envelope oscilloscope 
diagram is as shown in Fig. 2 and may be deduced 
from Fig. 1 as follows: between B and C only the 
direct sound is received by the microphone and at C 
the first echo arrives. Whether this arrives in phase 
or not depends on the distance BC and the fre- 
quency. If the distance BC is a whole number of 
wavelengths the echo arrives in phase and the en- 
velope is increased stepwise by an amount corres- 
ponding to the magnitude of the echo given in 
Fig. 1. If the distance BC is an odd multiple of the 
half wavelength the echo arrives in opposition and 
the envelope is decreased stepwise by the same 
amount. Thus from the echo time diagram is de- 
termined the relative phase of each echo at any given 
frequency and thus the diagram of Fig. 2 is con- 
structed. This diagram gives the transient at the 
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start of a tone and after a long enough interval the 
steady state response at the frequency in question. 


_ A similar procedure gives the switch off transient. 


Thus from the echo time diagram can be deduced 
all the essential information about the studio in- 
cluding the transients both starting and stopping 
and the steady state response. 
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Fig. la. Vector diagram. Relating Fig. 1 to Fig. 2. 
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Fig. la is the vector diagram providing the link 
between the echo time diagram (Fig. 1) and the 
envelope diagram (Fig. 2). The vector 01 represents 
the direct sound pressure. The vector 12 is drawn 
to represent the first echo at C both in magnitude 
and phase. The magnitude is taken direct from the 
echo time diagram and the phase is computed from 
the time delay BC between the direct sound and the 
echo allowing for the frequency of the signal. Simi- 


~ —— envelope 


Fig. 2. Envelope diagram. 


larly, the vector 23 represents the second echo and 
so on. The vector 02 represents the sum of the di- 
rect sound and the first echo and the vector 03 the 
sum of the direct sound and the first two echoes 
and so on. From this vector diagram is drawn the 
envelope diagram Fig. 2. 

Many specimens of experimental observations of 
both Fig. 1 and Fig. 2 have been taken and a de- 


scription of the apparatus and methods will be given 
elsewhere. If the auditorium is very irregular the 
echoes will be in random phase and no pattern will 
be detected in the echo time diagram, and the evalu- 
ation of the total response will have to be carried 
out by statistical methods [3]. There will be no 
favoured frequencies and the random variations 
will not be large. In the case of the rectangular 
room, however the characteristic frequencies and 
eigentones are evident as regularities of the echo 
time diagram. 


2.2.0. The rectangular room 


Suppose the dimensions of the rectangular room 
are «, B and y respectively, and suppose for simpli- 
city that the source is at one corner. Then the sound 
field in the room may be taken as that due to the 
source and its images in the walls of the room. Each 
image will itself have an image in the opposite wall 
and this again has its image and so on. The spacing 
of the images will be 2«, 28 and 2 in the respective 
directions and the whole constitutes a rectangular 
array filling all space. In Fig. 3 are drawn a few of 
the images in the horizontal plane. In a plane 2y 
vertically above this will be a similar set of images 
and this array is repeated at intervals 2y up and 
down. There is a great deal of regularity in this 
array and the regularity may be observed in two 
ways. The images are observed to be in straight 
lines and also in planes in various directions, 


J 


eB wn 
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Fig. 3. Images in horizontal plane. 


2.2.1. Line pattern 


The most obvious lines are those parallel to the 
walls of the room but it will be seen that there are 
also diagonal lines. Ifa straight line is drawn through 
any two images, an infinite number of other images 
also lies on this line. Take, for example, the source 
and the image n,2.« along, n.28 across, and n32 up. 
Then the image 2n,2« along, 2n,28 across and 
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2n,2y up are on the same straight line. The distance 
between successive images on this line is 


% = 2 (n, a)? + (ng 8)? + (ng y)?, 


where n,, n, and n, are whole numbers without com- 
mon factor. It is seen that on each line the images 
are equidistant and the distance is different for 
different lines. 

Now consider one such line and draw the echo 
time diagram for a microphone on the line ignoring 
for the time being all echoes from other images. The 
echoes will arrive at equal time intervals, the inter- 
val being x divided by the velocity of sound a, and 
the echoes will give sound pressures inversely pro- 
portional to distance (or less allowing for reflection 
loss and absorption). For steady frequency input, 
if the wavelength is x or a submultiple of x all the 
echoes will arrive in phase and add directly to a 
large value, and in the absence of absorption the 
value tends to infinity. For other frequencies the 
echoes arrive out of phase and mutually cancel. The 
pulse response can be evaluated as a Fourier series 
and as might be expected contains all frequencies 
whose wavelengths are x or submultiples of x. The 
components are also inversely proportional to time. 
These frequencies are not resonance frequencies in 
the usual sense although they are characteristic fre- 
quencies in that the response is greater than at non- 
characteristic frequencies. The frequencies are not 
eigentone frequencies. 

If the microphone is not on the line the echoes do 
not arrive from the first at equal time intervals and 
this deviation is more important for higher frequen- 
cies or shorter wavelengths. The later echoes do 
arrive more and more at equal intervals but as the 
echoes get smaller with time the first few are the 
more important and so the build up is much less if 
the microphone is not on the line, especially at 
higher frequencies. These line frequency modes are 
therefore dependent on position of microphone in 
a manner quite different from the eigentones proper. 


2.2.2. Images in planes 


Not only are the images in lines but also in planes. 
For example, if the source is on a floor corner there 
will be an infinite number of images in the ground 
plane as shown in Fig. 3. There will be a similar ar- 
ray of images on the horizontal plane 2y above 
ground and in the horizontal plane 4 above ground 
and so on. The distance between these planes is 2y 
and the number of images per unit area in each 
plane is 1/4«8. Similarly the planes parallel to the 
By wall of the room are distant 2 from one another 
and the number of images per unit area is 1/4By. 


The number of images per unit volume of space is 
1/8aBy, and clearly with any plane of images the 
number of images per unit area is proportional to 
the distance between planes. 

The oblique planes are not so easy to understand 
without the use of analytical geometry but one ob- 
lique plane which is the vertical plane through the 
floor diagonal is shown with parallel planes in Fig.4. 
The distance between the planes in this array is 


eae 
2/ V a: 45 Be ‘ 


Fig. 4. 


In Fig. 5 are drawn the traces of the vertical 
plane passing through the source and the image 
(3 X2a, 4 x26), and traces of parallel planes are also 
drawn. The distance between the planes in this array 


is Te ae 
4 \2 oe 
‘ yi) eae 
It can be shown by analytical geometry [4] that in 


the general case the distance between parallel planes 
is A=2/q 


2 2 
where g= (2) ea, +(e 
% B % 
n,, n. and ng being whole numbers without common 
factor. 
The number of images per unit area of such planes 


is \/8aBy—1/4a8yq and the mean distance bet- 
ween images in such a plane is 2 /aByq. 


. ° . . 
. . . . 
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| 2.2.2.1. Effect of images in a plane 


Now consider such a plane of images as heard at 
| a distance. If the distance is small only the images 
nearest the microphone affect the initial wave front, 
_ for the wave front from other images arrives too 
late. As the distance increases the effect of each 
_ image decreases but wave fronts from more images 
_ arrive nearly enough in synchronism. The effect is 
that the wave front from the plane of images be- 
comes flat and ceases to decrease with distance [5]. 

In Fig. 6, a, b, c etc. represent images in a plane nor- 
mal to the paper and circles have been drawn about 


a+ 
b+ 
es 

p Pp 
d+ q q 
e + 


Fig. 6. Wave front at various distances from a plane array of sources. Circles ‘‘p 
sections of the wave fronts from sources in the plane of the figure. Circles ‘‘q 
of wave fronts from sources above and below the plane of the figure. 


these points as centre to.represent the pulse waves 
by the time the wave fronts have travelled various 
distances. These circles represent the sections of the 
spherical wave fronts by the plane of the diagram, 
and other circles represent similar sections of the 
wave fronts from images not on the plane of the dia- 
gram. It will be seen that the total wave front from 
the plane of images rapidly becomes plane. The flat- 
ness of the wave front depends on the distance of 
the wave front from the plane of images compared 
with the distance between the images on the plane. 
It is to be noticed that although each image gives a 


—— pulse response 


A B 


time 


Fig. 7. Simplified equivalent echo time diagram for a system of planes. 


single impulse the plane of images gives a step func- 
tion, for in Fig. 6, behind the wave front the space 
becomes uniformly filled. The corresponding echo 
time diagram ultimately becomes like Fig. 7. In- 
stead of an impulse as in Fig. 1 at each echo there is 
a step corresponding to each plane of images, and 
the time between successive steps is the distance 
between successive parallel planes divided by the 
velocity of sound. 


After a long enough interval the room is filled 
with a series of plane wave fronts. The distance bet- 
ween successive planes is the wavelength of the 


6c 


are great circle 
‘q’’ are sections 
| 


lowest frequency which these wave fronts will sup- 
port. Since the distance between planes is 


In RPE 2 
A= 2/q, where q= y (2) bi (7) Es (28) t 


the lowest frequencies are 


a nm, \2 ty \? nN, \2 
VEG) 
where f is the frequency, 
a is the velocity of sound, 
71, MN, and n, are whole numbers without com- 
mon factor. 

Since all harmonics of these lowest 
modes are also possible, the frequen- 
cies supported by the regular arran- 
gement of the images in planes are 


PE ifaBs\ 25> (Bah Prge (tag\? 
cea badlern alten 
where k is any whole number giving 
the number of the harmonic of the 
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series of frequenciessupported by thesystem of planes 
in question. These are the eigentone frequencies [6]. 

Returning now to the echo time diagram, Fig. 8 
shows a simplified schematic illustration of the 
returns from one set of planes. At intervals corre- 
sponding to the distance between planes, there is a 
stepwise increase in the number of single image 
wave fronts arriving which ultimately gives a con- 


figuration like Fig. 7 


==ss oe 


——=- puise response 


i 


lig. 8. Simplified echo time diagram for a system of planes, 


From P to Q are the echoes from one plane. 
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soon sets of individual echoes appear in uniform 
time sequence. These give rise to transitory cha- 
racteristic frequencies which decrease inversely with 
time and soon vanish. At this stage, another kind of | 
regularity is observed among the echoes now very 
numerous. Large numbers of echoes arrive in groups 
very nearly simultaneously and the number in each 
group which arrives at nearly the same instant in- 
creases so as to compensate exactly for the smaller 
magnitude of individual echoes. This 
periodicity does not decrease except 
by absorption and gives the eigentone 
proper. 

Now consider the response at a fre- 
quency f = ka/A where as before 4 = 2/q 
and 


Vevey 


ill i 


From Q to R are the echoes from that plane plus the echoes 


from the next plane. 


At R begin the echoes from the third plane and so on. 


For a steady input those frequencies whose wave- 
length is a submultiple of the distance between 
planes build up without limit and, if the input is 
removed, continue indefinitely. These are the eigen- 
tones proper. 

If, now. reflection loss is taken into account, the 
situation is quite different. Consider how these ei- 
gentones are excited by the impulsive sound from 
the loudspeaker. Before the wave fronts can excite 
and maintain any frequency, the distances between 
successive wave fronts must be sufficiently equal in 
relation to the wavelength, or the wave fronts must 
be sufficiently flat. This requires that the distance 
from the images to the wave fronts must be large 
compared with the spacing of the images in the 
plane. In other words, time is required before the 
wave fronts are sufficiently flat to delineate higher 
frequencies and particularly so in planes which are 
close together and therefore contain few images per 
unit area. During this relatively long time, there will 
have been many reflections and consequently large 
absorption so that sound is attenuated before the 
higher frequency eigentones can build up. 

An approximate numerical figure can be obtained 
for the relative magnitude of the various eigentones 
as excited from a spherical source but before going 
on to this the total picture will be reviewed. Consider 
again the echo time diagram. Immediately after the 
input impulse the echoes arrive more or less at 
random depending on the relative positions of the 
source and the microphone. The echoes decrease in 
magnitude with time and increase in number and 


ris the distance between parallel planes 
of the group whose eigentone numbers 
are m1, mM, and ny. This group of parallel 
planes gives rise to frequencies whose wavelengths 
are submultiples of 4. In other words, this group 
of eigentones consists of a fundamental or lowest 
mode of wavelength % and frequency a/A and a 
series of harmonics of frequency 2a/A, 3a/A etc. One 
of these is the frequency f=ka/A=kagq/2. The 
distance between planes is 4 = 2/q and the number 
of images per unit area in each plane is 1/(4aByq). 
The mean distance between images on the plane 
is thus 2 «By q. 

Now consider the sound field at a point P distant 
s from one plane. In order to support exactly the 
required frequency the wave front from the plane 
of images should begin to arrive at P in time s/a. 
The wave front from the next plane should arrive 
d/a later and the next 24/a and so on. Actually the 
wave front will arrive exactly on time only if there 
is an image in the plane exactly opposite P, that 
is at the foot of the normal from P to the plane. As 
the mean distance between images in the plane is 
2y/ aByq the actual distance from the nearest image 
to the foot of the pormal from P to the plane will be 


The distance of this nearest image from P is 
) (s? + x?«8yq) instead of s and the deviation is 
Pe ip rere 
Taking a binomial approximation 
ds=x2 aByq/2s. 
Now the wavelength with which we are concerned 


is \/k = 2/kq and the deviation of the wave front at 
P as a fraction of the wavelength is 
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ORE 4s 
This is an angular deviation of 
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To obtain a reasonable approximation, suppose that 
0 is less than 7/3 or 60° so that the in phase contri- 
bution of the wave front is greater than 0.5 for the 
first few wave fronts. The average of x? is taken as 
1/5 so that 
¢ 2 ¢ 
2naPygq?k = 20 


2K 
SOs or Ss aii ° 
12s 6 2 


With this value of s the first wave front may be out 
of phase by 180° depending on the position of P 
relative to the plane. With oblique planes the images 
in one plane are not in alignment with the images 
in neighbouring planes so that if the error on the 
first plane is a maximum the error in the next planes 
will be less. 

In the case of oblique planes therefore the average 
of the first few planes gives an error less than 50% 
on the response of those first few planes. The planes 
still more remote will have progressively smaller er- 
rors. The responses of all the planes form a geome- 
tric series, the response of each successive plane 
being smaller in a constant ratio depending on the 
reflection loss. The sum of the series is proportional 
to the first term. In the case of the simple modes 
(1, 0, 0) (0, 1, 0) (0, 0, 1) the images in successive 
planes are in alignment and consequently the errors 
on all successive planes are cumulative so that the 
result is dependent on the position of P relative to 
the source. 

In the next section an approximation will be 
made of the different cigentones and here it will be 
sufficient to caleulate for a few particular cases the 
distance s and the corresponding time t= s/a that 
elapses before the eigentone begins to build up. Con- 
sider a room in which « =24m, 8=15m, y=7.5m. 
This represents a studio of 2700 m$ and the reverbe- 
ration time would be of the order 1.5 s. The results 
are tabulated in Table I for a few eigentones. 


Table I 
cio de | frequency distance s | time s/a 
| | 
10, (1, 0, 0) | 69 c/s 23.4m 0.071 s 
(10, 0, 1) 72 ¢/s 258 m 0.678 s 
(0, 1, 10) 220 c/s 


2400 m |’ 7.3 8 


Now consider a smaller studio say 6 4.53.6 m3. 
In this the reverberation time would be about 0.4 s 
and the results are tabulated in Table IT. 
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Table LIL 


eigentone 


eee hats frequency distance s time s/a 
| 

LOR (1, 050), 275 ¢e/s 13.5 m 0.041 s 

(10, 0, 1) 280 c/s 139 m 0.42 s 

(1, 0, 10) 460 c/s 378 m eh ers 


2.2.2.2. Discussion of numerical results 


In the larger studio in Table I the eigentone (0, 
1, 10) of frequency 220 c/s requires 7.3 s before the 
wave fronts are sufficiently plane to begin to esta- 
blish the eigentone pattern. By this time the attenu- 
ation by reflection is of the order 240 db. This eigen- 
tone is completely inaudible. In general terms it is 
clear that higher and oblique eigentones are heavily 
attenuated and only the lower frequency eigentones 
are of acoustic importance. 


2.2.3. Eigentone space pattern 


Before going on to state conclusions some state- 
ments on the spatial configurations of the eigentone 
pattern will be made. After a suitable lapse of time 
the studio will be filled with a large number of plane 
step wave fronts. These wave fronts travel norma] 
to themselves at the speed of sound. Corresponding 
to each value of 


VEY +GY+EY 


there will be in all eight sets of wave fronts. In each 
direction of four directions there will be two wave 
fronts travelling in opposite directions, thus creat- 
ing a standing wave pattern. The total pattern 
consists in essence of a set of parallel travelling wave 
fronts at equal intervals plus a similar set travelling 
in the opposite direction plus all the reflections of 
these in the walls of the room [7]. These together 
constitute the spatial pattern corresponding to that 
eigentone, and the pattern is complicated for large 
values of q. 


2.2.4. Excitation of individual eigentones 

If it were required to excite immediately the 
simpler eigentones, for example, the (1, 0, 0) mode, 
it would be sufficient to drive one end wall of the 
room parallel to itself at the appropriate frequency 
but to excite a mode with high value of g would 
require a very complicated system of driving planes. 
Failing the simultaneous drive of all the appropriate 
planes time is required to establish the pattern and 
as is seen from the preceding, a relatively long time 
for high values of q. 

Another way of looking at this is to consider the 
spherical impulse source not only as containing all 
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frequencies equally but also as containing all orien- 
tations equally. It is so to say impulsive not only in 
time but also in space. Since each eigentone has not 
only a specific frequency but also a specific orien- 
tation it follows that the degree of excitation of any 
eigentone depends not only on the frequency spec- 
trum of the source but also on the polar diagram 
of the source. 

Practical sound sources are of course not spherical 
but decidedly directional and the relative excita- 
tion of the eigentones depends largely on the orien- 
tation of the source or sources. Nevertheless any 
practical source, however directional, sounds like a 
spherical source at large distances and highly ob- 
lique high order eigentones will be only faintly 
excited even with highly directional sources. 


2.3.0. Conclusions 


We are now able to make a general picture of 
what happens in a rectangular studio when pulse 
excited and to deduce from this in general terms 
what happens when steady tone is switched on at 
any frequency. The results are best considered in 
three stages. 


2.3.1. First stage 


In this stage, the echoes are very few and large 
and there is as yet no periodicity imposed by the 
room. The interference patterns produced depend on 
the relative positions of the source and microphone. 


2.3.2. Second stage 


In this stage, the echoes are more numerous and 
smaller and the line pattern is noticeable giving rise 
to frequencies which however are not eigentone fre- 
quencies. These decay inversely as time whether or 
not there is absorption. These frequencies are in 
essence those causing the bulk of the audible effect 
of the studio. The rate of decay is characteristically 
not exponential. At the end of this stage the sound 
has generally decreased by thirty or forty decibels 
from the peak values. 


2.3.3. Third stage 


In this stage the echoes generally are arriving in 
groups nearly in synchronism from the various 
planes of images. These establish the eigentone 
pattern proper which is similar all over the room in 
distinction from the previous stages where the 
patterns were largely local. Owing to the time taken 
to establish these patterns, in a studio with usual 
absorption, the higher frequency eigentones are low 
in level and the highly oblique modes are virtually 
absent. The die away of this stage is characteristi- 
cally exponential. 


Steady-state characteristics 


A true resonance can only be built up by a large 
number of echoes. This is implicit in the basic idea 
of frequency. As stated for example by Gaxpor [8] 
if frequency is known precisely epoch is completely — 
indeterminate and conversely if epoch is known pre- 
cisely frequency is indeterminate. The impulse is 
precisely known as-to epoch and contains no specific 
frequency. — 

It requires an infinite series of impulses to de- 
termine a frequency precisely and more generally a 
sharply defined resonance requires a large number 
of echoes in regular sequence to define it. There is 
however another phenomenon which resembles a 
resonance in some respects but which is in effect 
quite different. With any particular location of 
source and microphone it may happen that the first 
few larger echoes give a zero or nearly zero total 
response. If the steady state characteristic is on a 
decibel basis, as is usual, at such frequencies there 
will be a deep narrow crevasse on the curve. This 


Fig. 9. 


crevasse is not a resonance phenomenon in the usual 
meaning of the term. In Fig. 9 is shown a small sec- 
tion of a steady state characteristic for a typical 
studio. At P is a crevasse of this type. These cre- 
vasses are a prominent feature of steady state cha- 
racteristic curves and it is on the other hand diffi- 
cult to find any true resonance feature on the curve. 


3. Mathematical section 
3.0. General 


Suppose at the point P a velocity source of such 
a nature as to give at unit radius a positive unit 
pressure impulse [2]. This will be taken as unit 
source. From the source a pressure wave travels 
radially outwards with velocity a. The pressure at 
radius r will be 


P(t) =+ 3(:-1) (1) 


where 8(t—r/a) is zero for all values of t except r/a 


and such that 
8 ( : ja = 2 
ari t - ( ) 
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In the notation of van per Pot if V(p)= P(t)*) or 
V(p) is the “image” of P(t) 


Spe iP 1 — pri/a 
Vip)=p |e _ 3(: ze —pe *. 
é (3) 
Since the input is represented by p this may be 
taken as 
V(p)=pZ(p) 
where Z(p) is a transfer impedance ratio and the 
response to any other input 
Pol) Voip) is, 
V(p)= Vo (p)Z= Vo(p) ——e (4) 
It is easily seen that this is the “image” of a solu- 
tion of the equation for a velocity potential 
do 
de 
Since pressure is a scalar quantity 


—pria 


— a*7V20=0. 


V.(p) = > Vo(p) = ar Pele ov. %() 
gives the solution for a series of sources. 


3.1.0. The rectangular box 
3.1.1. Formal equations 

Consider a rectangular box whose edges which 
meet at a point are taken as axes of rectangular co- 
ordinates and are «, B, y respectively. Take as source 
the point P (x,, y,, z,) and observe the resultant field 
at Q (%o, Yo, 2). The field will be taken as due to the 
source plus all its images in the walls of the box. 
Thus for unit impulse source 


Z(p) ay a (8) 
For any input P,(t) = Vo(p) the response is 
V(p) = Vo(p)4 (9) 


and if P(t)= V(p) the response in terms of time 
is formally 


I Vo(P) 

ist ba by teas pt _O\F 7) 
P(t) = V(p) gay fe é Zdp _ (10) 

where the integral is the Bromwich integral. 
The response to unit impulse is 
1 
ee eS pt 

P(t) Fare Zdp. (11) 


eC 
3.1.2. Solution of the formal equation 


These equations (10) and (11) may be solved in 
various ways. The expansion in terms of the singular 
points of Z(p) corresponding to Heaviside’s partial 
fraction expansion gives results in terms of ex- 
ponentials. The imaginary exponents of these terms 
correspond to the critical frequencies implied in Z(p) 
and thus imposed by the rectangular box. The ex- 
pansion thus consists in finding the singular points 


of Z(p). Now considered as in circuit technique 


Z(p)= — exp [—(p+v)r/a] consists of a sum of 


time vectors of various phases. If an ordered set of 
these vectors coincide in phase so as to give an in- 
finitive sum, this constitutes a singularity of Z(p). 
Now the source and its images constitute a rect- 
angular array filling all space and these may be 

grouped for summa- 


2| tion either in lines or 


mn 
B "| in planes. As stated 


Pacxp 2 Vey tay Arya)? + (Yo +1 + 278)? + (% + % 
V(p) = sf 
® V (% ay t 270)? + (Yo t ¥y + 27g8)? + (% + % + 275 y)? 


— » P- exp (—pr/a) (6) 
- 


Z(p) =y'ee Cah) (7) 


r 
This is for the case of zero reflection loss or absorp- 
tion. If reflection loss is taken as e—? per reflection 
then a factor exp [— (r, +1r,+173)e] is required. 

If an equivalent loss as e—¥! is taken to secure 
approximate results then 


and 


*) Hither V(p) = U(t) or U(t) = V(p) is an abbreviation for 


Vip) = p{ePuwat 
0 


which is well known to be equivalent to 
1 V (p) 
= pt 
U(t) al “ dp. 


Cc 
C being the Bromwich contour. 


before the grouping 
in planes gives the 
eigentone solutions and the grouping in lines gives the 
expansion in terms of the Fourier series for the sets 
of echoes. Having found the singular points the re- 
mainder of the formal solution consists in evaluating 
the residues corresponding to equations (10) and (11). 


3.1.3. Symmetry 

The rectangular array of images is symmetrical 
about the six planes which are the walls of the box. 
The development of this symmetry meets the 
condition of zero velocity normal to the walls of 
the box. Immediately after the initial pulse, how- 
ever, the sound field in the box shows no symmetry 
or order. It is only after several reflections have 
taken place that order and symmetry appear and 
consequently specific characteristic periodicities. 
By this time x,, x, etc. are small compared with 
2 V/(r,«)? + (728)? + (r3y)2=m so that Z(p) can be 
expanded by the binomial theorem approximately. 
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Taking 2r,«/m, 2r,8/m, 2rzy/m=L, l,, 1, respec- 
tively as the direction cosines of the line from the 
origin to the image (217, «, 27,8, 2r3y) 


Z(p) a exp[—(p+u)m/a] | gy 9 


m 
where 
l l, x L 
Q=cosh es 1 P alee oon eae!) 21 
ip y ee 
cosh wale 2Y2 oh (P ss 3 16 eee) 3% 


neglecting terms in 1/m?. 

In the final development (p+) takes some 
imaginary values so that these terms become the 
well known cosine space distribution terms of the 
field. To avoid carrying these terms throughout, 
make x, =*.=Y¥,= Y2=2,=%,=0. This corresponds 
to placing the source at one corner of the box and 
listening at the same corner. In this way all the co- 
sine terms in Q are unity or in other words all pos- 
sible frequencies can be excited and all excited fre- 
quencies can be heard. 

With source or listening point anywhere else some 
component frequencies will not be excited and some 
excited frequencies will not be heard. Making this 
simplification Z(p) becomes (12) 

6g OP le tw)laly 2ra)*+ 2rg8)+ 2rsy)*} 


Z(p)= ni ee = 
V@r, 2)*+ (27,8)*+ (2157)? 

In this case the source and its first seven images 
which are coincident with it together form a virtual 
source eight times greater and similarly the micro- 
phone is multiplied eight times. This accounts for 
the factor 64. 

It is to be remembered that this symmetry is not 
evident during the first phase immediately after the 
input and this will be investigated later. 


3.1.4. Expansion in terms of line frequen- 
cies 
Taking eq. (12) and grouping in lines 


Z(p) = >= [—(p +u)2 km/a] 


64 2km 
m k 


where m? = (r,«)? + (r,8)? + (rg)? and r,, rj, Tr; take 
all sets of positive integral values without common 
factor, including zero for any two in any set. For 
example r, =r; =0, r,; =1 gives the term 


exp [—(p+u) 2ka/a] 
Qha 
Now 


k 
d Z 
oe a - ax [—(p +u)2km/a] 
m k 
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(13) 


P L 32 Oy 
Wa 1—exp[—2 (p+ u)m/a] 


m 


yn — ut , kant 

hte e : 

: Tie (i+: cos — — J. 
: k 


me 


Now if f(p)=h() =p | eP*h(@) ds 
0 
d_ f(p) pri ie 
—pa=——— =p |e Pit k@ias 
dp BZ | 
=t h(t) 
eve kant 
Thus vipi=os > iia (i+) om a 


where m=V/(r,a)* + (r28)?+ (rsy* (14) 
and r,, T2, T3 are integers without common factor. 

This we shall call the Fourier series expansion of 
the response in terms of successive echoes. It is seen 
that the factor 1/t makes these periodicities fugitive. 
The frequencies are not eigentone frequencies. 

It is now to be observed that as heard from the 
origin the echoes in all these series arrive at equal 
time intervals so that the result of eq. (14) is cor- 
rect. At other points where in general the listening 
point is not in the line of images making up the 
series, the echoes will not arrive from the first at 
equal time intervals but the deviation from equality 
rapidly decreases with time. The deviation is more 
important at higher frequencies where a given time 
deviation represents a greater phase deviation. Thus 
a longer time must pass before the higher frequencies 
develop, and owing to the factors 1-1 e—¥* the 
higher frequencies are reduced. At specific locations 
for source and listening point, where the listening 
point is on the line of images concerned, this re- 
duction no longer holds. 

In general these line frequencies are dependent 
on the relative positions of source and listening point 
in a manner quite different from the eigentones 
proper. The eigentone field apart from the arrange- 
ment of nodal and loop planes is uniform over the 
room whereas the line frequency field is decidedly 
a function of position of source and microphone. 


3.1.4.1.Closer approximation of absorp- 
tion 

In the case of the line periodicities a closer ap- 
proximation may be made for the absorption than 
the assumption of attenuation as e—#*, Suppose at 
each reflection a reflection loss as e— °. Then clearly 
successive terms in the line series will be attenuated 
as exp [—(n, +n, + n5)¢]. Thus from eq. (12) 
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eD Dy exp (—2pkm/ja) 
64 P 2 km 
m k 


-exp[—(n, +n, + 15) pk]. 
Then finally (15) 


>; — ; : k 
P(t) = 64 25k: Beas mt) ( + 2 cos aol 
m k 


This is the same as (13) except that u is replaced by 
nae (mtmtns)ap 
ETO i oy 2 2 2° 
(m, «)? + (m1.8)? + (73 y) 
3.1.4. 2. Sinusoidal input 
The response to sinusoidal input could be deduced 


(16) 


. * ct 2) 
from eq. (12) by putting u)=sin wt = 10 52 + we 
so that 2 

@ 
AV Wesecimmpe Ata} 


It is, however, easier and more illuminating to 
proceed as follows. From eq. (6) can be deduced the 
magnitude and epoch of echo. Suppose the input is 
sin wt and the frequency corresponds to one of the 
line frequencies. Then depending on the positions 
of the source and the microphone, certain of the 
echoes will arrive at equal time intervals correspond- 


ing to 
2m 2 
“aa V (4)? + (8)? + (gy)? 


The effect of these will be summed since they arrive 
in phase. The remaining echoes will arrive at random 
so far as this frequency is concerned and their effect 
can be estimated statistically as will be seen later. 

In the simpler case of the source at one corner 
and the listening point at the same corner and 
omitting the direct sound, the series to be added is 


k=o 
) 1 1 
k=1 


There are, however, in this case eight sets of 
images in lines which also contribute, the co-ordi- 
nates of these being +2n,a, +2n,8, +2n,y, thus 

Q=—(4/m).log (l—e—"P). (17a) 

The next step is to obtain an approximation for 
the sum of the random phase echoes for the sake 
of comparison with these values and later with the 
eigentone values. This will be pursued in the next 


paragraph. 
3.1.5. Random phase summation* 


Suppose the phase of the various echoes at any 
given frequency is taken as random. Then the sum 


(17) 


* The argument which follows is derived from HORNER’S 
paper (reference [3]). 


of the components in phase with and in quadrature 


‘with a given arbitrary vector are 


x= a, cos 8, + a, cos 9, + ..., 
¥y =a, sin 6,+ a, sin 0+ .... 
therefore x? + y?= Za,?+ La, a; cos 9s cos 
+ Da, a; sin 6, sin 0;. 
The mean square value is the mean of «2+ y? over 
all values of 6,, 9,,... 9, which is 
27 27 an 
( 


1 2 
way [Af a... [et +72) ae, = as 
6 0 


fy 
since the integration of all the cosine and sine terms 
gives zero. 

Thus the vectors in random phase add as if in 
quadrature to give the mean square value. The 
R. M. 8. value is the square root of this. This is on 
the basis that the numerical value of the echo is in- 
dependent of the phase. This is not strictly true, 
especially at low frequencies, but the result is prob- 
ably reasonably accurate since although the mag- 
nitude and phase are related they are not related 
in any regular way as for example if increase of r 
always corresponded to an increase of 6,. 

Now a; may be expressed as 

as oe e trs/a, 
Ts 
Here p is an average decrement figure related to the 
average reverberation time of the studio in question 
by the equation p= 6.9/r where + is the rever- 
beration time. 


In the case of a rectangular studio, two cases are 
to be considered. 

1. When the source or the microphone is near a 
corner of the room. “‘Near’’ means that the dis- 
tance from the corner is small compared with the 
wavelength considered. 

2. Other cases. 

Case 1. In this case the source and the first seven 
echoes are close together and make virtually a single 
composite source. The room corner acts as a horn. 
In this case the number of images per unit volume 
is n= 1/8 V= 1/8a8y where Vis the room volume. 

In time t the sound has travelled a distance at =r 
and all images inside a sphere of radius r and vo- 
lume $xr® will have been heard. The rate of arrival 
of echoes is 


dee pe 
Aravehekeic Shs 4 


per unit of length. 
Hence the mean square sum of all the echoes is 


xe ih gape PHr lady = Se Parla, 
ue 


mr? mr? 


Ma 


(18) 


Mm 
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Taking the direct sound at unit distance as unity 
and including in this direct sound all the group of 
the composite source, “‘r,”’ is the range of the nearest 
other echo. Since pu is small in practice r, need not 
be accurately estimated. 


Case 2. In this case 


9 
ge Soren e— BN /a 


PL Te 


again taking direct sound at unit distance as unity. 


(19) 


3.1.6. Expansion in terms of eigentones 


In the case of the source at the origin the images 
will be at the points +2r,«, +2r.8, + 2rey. 

The number of images per unit volume is 1/8 «8 y. 

Now consider the plane 


TH, May. Mek 

o e x 
where 71,, n,, nz are integers without common factor, 
It is easy to see that an infinite number of images 
lie on this plane. 


Let 


o-(2y+(3)e(3) 


then I, m, n, the direction cosines of the plane are 
n,/«q, M/Bq, Ns/~q respectively. 

Now consider the point r,2«, 7228, 732 y. 

The distance of this point from the plane is 


s=2Ir,o+ 2 mr,6+2nrgy= — (1, +12Mg +13 Ng). 


Since r, 73, Tz are any positive or negative integers 
1,1, + T.N, +733 can have any integral value and 
the lowest value of s is thus 2/q. 


Through this nearest point to the plane the equa- 

tion of the parallel plane is 
Ix -+my+nz= 2/q. 

This plane also contains an infinite number of 
images. Thus the images may be grouped in a sy- 
stem of parallel planes, distant 2/q from one another 
as expressed by the equation 
Ix-+my-+nz= 2k/q where k is any positive integer. 
ae uM Da a soltad 2k 

BF bol BG orch Mt bong 
with the same value of q and therefore the same 
spacing, there are four sets of such planes 


TX | MY , Mgt _ 2k 
aq ~ Bq 19. 4" 
These are actually one plane plus all its reflections 
in the walls of the box. 
There will be another similar set of planes for each 
value of q, that is for all sets of values of n, ng, n3. 
The distance between planes is s=2/q and the 
number of .. ages per unit volume is 1/8aBy and 


(20) 


thus the number of images per unit area of each 
plane is 

n=1/4aByq. (21) 
The average distance between images in the plane 
is 2//aByq. 

In general the images in a plane will lie at the 
corners of a parallelogram lattice. The normal to the 
plane drawn through one image does not, in general, 
pass through images in other planes. That is, the 
images are not, in general, aligned normal to the 
planes. 

The maximum distance of any point in a plane 
from the nearest image is taken as 


1 =Vabyg. (22) 
It is now required to find an expression for the field 


due to the images in one such plane. 


3.1.6.1. Field due to images in a plane 


Consider first a plane with the discrete images 
replaced by a large number of small images but with 
the same density per unit area. Then by direct inte- 
gration the pressure at distance s normal to the 
plane is 


oxp(—2 yes) oa, 
ft Veta. “4aBbyq 


V(p)s= | p 


ait 
2Qare Psa 


This represents a plane wave of magnitude indepen- 
dent of s. 
If instead of impulsive input the input is 
2 


ye 
cos ates +o 


then the response is 


? P 
V(p)z= Pp’ devia (—2yr+s) 2nrdr 
Pp? +o? Vr2+ s? 4aByq 
— ps/a 
2 are chia F 1 f P oul (24) 
4eByq © p? +o? 
_ 2an 1 : s H 
Tidapye = sino (-— +). (25) 
If the limits are taken r, to oo 
2an P 
VO ag oP | “ V8+ ri). (26) 


In front of a hole in the infinite plane the response 
remains equal in magnitude but is delayed a distance 


is = Vs?+17r2—s. 
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3.1.6.2. Approximation for discrete images 


_ Now consider the error involved in taking this 
integral for the uniform plane source instead of the 
correct sum for the discrete images. 
~ For large values of s the difference is that behind 
the wave front the pressure is uniform in the inte- 
gral case whereas in the case of the discrete images 
it is non-uniform. 
_ Itis, however, precisely this non-uniformity which 
constitutes the recurring pattern giving rise to eigen- 
tones other than the one now considered. 
For small values of s the wave front is no longer 
exactly plane and is delayed except exactly oppo- 
site an image. 


3.1.6.3. Case of no absorption 


In the case of zero absorption s can be made 
arbitrarily large since the field from the plane does 
not diminish with distance. Ultimately the box is 
filled with a whole series of plane wave fronts and 
these constitute the eigentone pattern. 


3.1.6.4.Time taken for eigentone to de- 
velop 

As stated above, for small values of s the wave 
front is no longer plane and the importance of the 
deviation from flatness depends on the wavelength 
of the sound being considered. When the wave front 
is sufficiently flat, successive wave fronts arrive 
sufficiently in phase to add and give a large cumu- 
lative result. An approximation will now be made 
to the time required for the wave fronts to be suf- 
ficiently flat to support.any specific frequency. 

From any point P distant s from one plane, draw 
the normal to the plane. The distance of the nearest 
image in the plane from the foot of this normal 
cannot exceed VaByq by (22). In the general case 
the distance in the next plane will differ and the 
distance can be written 


rT, =xy aByq where x<1l. 
Hence the distance from P to the nearest image 


Approximating 

Raber as 
2s 

Now for any value of q there are, as stated before, 
four sets of planes of images. From each plane of 
images there is a wave front in both directions so 
that there are in all eight sets of wave fronts. In each 
direction the two wave fronts travelling in opposite 
directions set up a standing wave pattern. This con- 
stitutes the eigentone pattern. The first wave front 


(27) 


to be considered as sufficiently plane to build up the 
pattern will consist of the sum of two oppositely 
travelling wave fronts and so on. The average value 
of x? may be taken as 1/3 so that 
aby 
Gs 
Now consider a frequency in the harmonic series of 
eigentones for which the value q applies. 

The wavelength of the fundamental is 2/q and 
the wavelength of the kth harmonic is thus 2/qk. In 
order that the wavefronts may build up at this fre- 
quency and wavelength, 5s must be small compared 
with 2/qk. If, for example, 8s is 1/6 wavelength on 
the average, then for the first few planes the actual 
phase error will vary between 180° and zero so that 
for q large and consequently successive planes near 
together the sum of the first few terms will tend to 
have a zero in phase component. It is reasonably 
certain, therefore, that this value of Ss is near the 
upper limit. 9 

Taking this value of $s < 6qk 


Reeng 24 ane Be 
6s 64Gb and s > nt: 


NI 
SSaver, == 


(28) 


(29) 


This value of s as stated above is near the lower 


limit for s. 

If 8s is taken as 1/12 wavelength the error varies 
from 0 to 90° averaging 30° with an inphase defect 
of about 13%. 

Thus s may be taken greater than «af yq?k/2 
and less than «8 qk. . 

From this the time that must elapse before an 
eigentone begins to build up is greater than 

aByqrk 


a (80) 


3.1.6.5. Summation of eigentone series 


From eq. (23) the pressure from the plane is 


The next plane is more distant by a distance 2/q 
and so on. Allowing for the series of wave fronts 
travelling in the opposite direction, the sum is 


On Pa (1 40 2PM + 9 APlary....) 
“By 


The time function corresponding to this is 


1 ar / 1 
eee noe te pia) ee 8 dp, 
es gail cie3 e p (l—e—P/24) ‘ 
z 
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The zero frequency term tends to infinity owing 
to the nature of the impulse input and the total 
absence of absorption, and the response at other 
frequencies is 


sinkagr s s 
P=» aBygk (— a) a 


k 


That is, the fundamental of each set of eigentones 
corresponding to any value of q is inversely as fre- 
quency and in each harmonic set the components 
are also inversely as frequency. 


3. 1.6.5.1. Summation 
sorption 


allowing for ab- 


In the first instance allow absorption as e— 4 
where pu. is a general average figure which may be 
derived from the average reverberation time of the 
studio as follows: The reverberation time is the time 
taken for the sound pressure to fall 1000 times so 
that if + is the reverberation time 


1 6.9 
UT — e— 6.9 Sages 
eu T0007 © and uy = mace 
Then eq. (31) becomes 
ax 3 1 
Se or em 
he) aByq ° 5 1—e— (2p + w)/aq 
ka 
1 Siivelal sae einen s 
and Poy x He i ayo 
I 
approximately. 
2h 
Nowe == by (80) so that the numerical 


magnitude of the term of frequency kaq/2 is 
eH aBYPk/24 


Cay ee 2 
aByqk ; Fe 


3.1.6.5.2.Second approximation for re- 
flection loss 


In the plane 


MX | MY , M32 2n 


= + 

% 8 Y q 

the important contributions to the wave front in the 
box are from the images angularly near the foot of 
the normal from the origin to the plane. The equa- 


tions to this normal are 


Hence r, + r, +r, is approximately 
Cy LL | n Nn. 
ntntn=90 (gt Be +7) 


(where s is the normal distance of the plane) 


where 


n, \? Thy \ 2 ny \2 

ad gt=(H) + Gee. 

The reflection loss is taken as e° per reflection 
and for r, + r, + r,reflections 

e Titre thse — 9—aget/2q_ .—yt 

In equation (32) replace u by ag e/2q. 

Thus the numerical magnitude of the term of 
frequency kaq/2 becomes 


—“Byqkge/t os 
Soe Oe aa (33) 
aBygk 
3.1.6.6. Response to sustained sinusoidal 


input 

This could be calculated from the general equa- 
tions but it is easier to proceed as follows. Suppose 
the input is cos wt starting abruptly at zero time. 
Suppose o = kagqn, that is, the frequency of the kth 
harmonic of the series of which q determines the 
system of planes. 

The response from a plane array of images is 


2an osha / 1 Poe 
SS SS PS. a ¢———— 
LO) ead (+> aap 2 from (24). 
Allowing for absorption as e— #4 
V (page. GP he 


4aByq p?+o® (p+) 
Up to the point where s > «8 yq?k/2a the echoes 


add at random and after this add in order to give 
the total 


an Pp? 1 
abyg p?t+o? pty 
(Ite 2(P+eN/a9 4 o— 8p e ag ee 


on FP 1 pte 
~ aBygq p?+o% pty 1—e—2(p+H)/aq" 


“ oP Tr u.) 8/a F 


The steady-state response from this is 
an e Hsla 


= aya | ———2Heg 


uCOS wot—w sinewt 


w? + p? — + (34) 


P(t) 


Taking » small compared with this is appro- 
ximately 


an 6 PAs sin kagqr 


, etn (35) 
aByq y¢—2u/aq kagz 


P(t)= 


1 e He/a 
ae 2) 1—e—2H/aq 


sin kagqr 
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and inserting values of ». and s 


Sige Le 8g? sin kaqr. 


(36) 


3.1.7. Summary 


The results for the simple case of the source at a 
corner and the listening point at the same corner 
may now be summarised. 

There are three results and in each the response 
is given as a fraction of the direct pressure at unit 
distance from the source. The direct sound pressure 
is taken as the sound from the virtual source eight 
times greater than the actual physical source owing 
to the seven coincident images. 


3.1.7.1. R. M.S. value on basis of random 
phase 


Rit) 


4Vu 
where r, is the minimum radius outside which the 
summation is taken. 


fe er (18) 


3.1.7.2. Line frequencies 
Lit a log (I—e_”°) by(17a) where 
m? = (n, x)? + (m8)? + (gy)? and n =n, + ny + Ny. 
The frequency is ka/m. 


3.1.7.3. Eigentone frequencies 
I e— Vkagel4 


EQ) = Vek 1.80? by (36) 
where V= room volume =«8y 
gore fein Eos Nitty \* Mz \2 = My | My | Ng 
eo (2) ae) 07 RCI are Sina rn 


r is the order of harmonic in the series based on q so 
that the eigentone frequency is kaq/2. 


3.1.8. Numerical results 


To see clearly the significance of these results, 
the results are computed for two average studios 
at various frequencies. These are tabulated in 
Table IIT. 


3.1.9. Effect of position of source and 
microphone 


The effect of the positions of source and listening 
point differs in the three phases considered and can 
best be taken into account by a general review of 
the successive phases of the sound field after a signal 
is sent from the source. 


3.1.9.1. First phase 


During the phase immediately after the signal 
starts, there is no order or symmetry about the field 


WAVE PATTERNS IN STUDIO ACOUSTICS 


Table III 
studio | a{m]/B(m]/y{m]| t[s] vis] ¢ 
eae | fl | 
studio A | 24 | 15 | 7.5 | 1.5 | 46 | 0.12 
studio B | 6 3.8} 2.4 1 0.4 |.17.25; 0.13 
Higentones 
studio | n | r | | | frequenc | ree 
Hei ie 2 Ns | | a yI pressure 
Be Way | ole o | 20 | 1387.6 | 0.025 
eases OL ste ol Ths STB WAGS™ A) ORS Sai 
A ON LORR pe 20 440 Hr OF59SalO=—* 
A Liman: 20k omailaul 4. £40 0 
B 1 Oar) 10 244 0.015 
B LOMA OM 1 283 0.024. 
B 0 0 1 10 688 0.007 
B Desi Olmert pkO 1 688 6.9 x 10-8 
Line frequency values 
sates | | relativ: 
studio | ny | No | N3 k eee ees 
Ay cha |o! o |20:! 137.8 | 0.0136 
A J 1 Leal) 127 0.0069 
B 1 Onn 20r 4 10 274 ' 0.0525 
B 1 1 Lele 220 0.0217 


R.M.S. values (excluding direct sound) 


studio A 
studio B 


0.0012 
0.0195 


and the successive echoes add or mutually cancel 
in a manner highly dependent on the relative posi- 
tions of the source and microphone. This gives rise 
to interference patterns showing apparent resonan- 
ces which, however, are a function of source and 
microphone position relative to the room rather 
than a function of the room itself. A little later the 
R. M.S. random sum of echoes tends towards a 
value less dependent on frequency but the sum 
of this with the direct sound still remains highly 
variable with frequency and relative positions of 
source and microphone. 


3.1.9.2. Second phase 


In this phase the line frequencies begin to be evi- 
dent in certain relative positions of source and mi- 
crophone. The dependence on position is much less 
than in phase 1. These frequencies are related to the 
room dimensions but are not eigentones. They de- 
crease inversely as time and not exponentially as do 
the eigentones proper. 


3.1.9.3. Third phase 


In this phase the eigentones begin. These are 
dependent on position in the sense that there is a 
regular pattern of nodal and loop planes. They are 
independent of position, however, in the sense that 
any nodal point is as good as any other for the 
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microphone position. The origin is a convenient ana- 
lytical position because it is a nodal point for all 
modes. 

As remarked earlier, the position of the source has 
another effect. If, for example, the source is at the 
geometric centre of the room then the number of 
images is eight times greater and the planes are at 
half the distance from one another. Thus, for a fre- 
quency kaq/2 with the same values of n,, n, and n, 
as before, the number of images per plane is four 
times greater and the spacing half so that for the 
same frequency the order of harmonic is halved. 
The fundamental and odd harmonies of the original 
eigentones are not excited. The eigentones are as if 
for a room of half the linear dimensions of the actual 
room. The distance s is thus halved and the attenu- 
ation of the higher eigentones is less. 

Space does not permit a more detailed considera- 
tion of the effect in general of source and listening 
positions but the methods outlined above can be 
used to compute the effect with any specific ar- 
rangements. 

Appendix 
Velocity source required for a given pres- 
sure wave 


The equation to be solved is 


d?® 
ae 7 eae Oe (1) 
From spherical symmetry © is a function of r only 
so that the equation may be written 
d? d? 
Aa ee a =O (2) 
Following the method of van per Pot, multiply 
by e—P! and integrate between limits 0 and o. 


Denote p [oP (®) dt by (p) so that Y(p)= ® (t). 
0 


Then integrating by parts 


2 . 2 d? » d 2 
p?(r'¥)—a@ a3 (r¥)= Pai (rD)+ p?(r®) | t=0. 


Now initially all is at rest so that . (r}) =O 1 D0 


and 


‘ ad? 
and 
r¥=A.e 1/4 4B. ePrla (3) 


where A and B are independent of r. 

Now suppose at radius « a spherical shell con- 
strained to move with radial velocity U)(t) = t)(p) 
the displacement supposed small. This shell re- 
presents the velocity source. Since the shell is con- 
strained what happens inside the shell has no effect 


outside the shell and consequently only the out- 
going wave from the shell need be considered. 
Thus r¥= Ae Prla (4) 
where A is to be determined. 


Now if P(t)=- V(p) represent the pressure and — 
U (t) ==t(p) represent the velocity 


ip=Geo—(a+E)4 or ©) 


r2 ar 


A 


1 be 
Now Pale t ZA pa/a 


and this is equal to i) (p). Now taking zero of time 
not on the shell but from the time at which the 
wave would have passed the centre of the sphere 
the initial velocity is 


‘(ple wt 


ip) =—(9+ PY 4. (7) 


so that 
Also suppose 


1 : 1 


a 
Then A=— Vo(P) 
Pe 
a 
and is(?)= (Sat =) off) (8) 


t 
1 1 
and U0) — Sas P, (t) + =a Pol dt. (9) 
0 


This gives the velocity input at radius « to generate 
any particular pressure wave. 
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PROBLEMES DE RAYONNEMENT EN ACOUSTIQUE DU BATIMENT 


par M. JACQUES BRILLOUIN 


Directeur du Centre d’Acoustique Technique, Paris 


Sommaire 


On montre, par des exemples tirés de l’expérience et par des considérations théoriques que 
Yamplitude moyenne des vibrations des parois d’une salle ne peut renseigner sur leur rayon- 
nement, et par suite sur la voie empruntée par le bruit pour pénétrer dans la salle. Une paroi par- 
courue par des ondes de flexion ne rayonne de l’énergie que si la vitesse de propagation de ces 
ondes est supérieure & celle du son dans l’air; son <indice de rayonnement> peut étre trés dif- 
férent selon qu’il s’agit de vibrations forcées ou de vibrations libres. On esquisse une théorie du 
plafond suspendu qui rend compte de certaines observations. Enfin on montre que les phéno- 
meénes transitoires présentent des aspects différents de ceux du régime permanent, et on expose 
les résultats de quelques expériences relatives & cette question. 


Zusammenfassung 


Erfahrung und theoretische Uberlegungen zeigen, dafs die Kenntnis der mittleren Schwin- 
gungsamplitude einer Trennwand nicht ausreicht, um die Schallabstrahlung in einen Raum zu 
berechnen. Es 1a8t sich daher auch nicht voraussagen, auf welchen Wegen Stérgerausche in den 
Raum gelangen. Biegewellen in einer Trennwand strahlen namlich nur dann Schallenergie ab, 
wenn ihre Ausbreitungsgeschwindigkeit gréBer ist als die Schallgeschwindigkeit in Luft. Das 
,, Strahlungsvermégen‘‘ einer Trennwand kann daher sehr verschieden sein, je nachdem ob es 
sich um erzwungene oder freie Schwingungen handelt. Es wird eine Theorie der hingenden Decke 
gegeben, die einige Beobachtungen erklaren kann. SchlieBlich wird gezeigt, da8 Einschwingvor- 
giinge Besonderheiten zeigen, die man im stationaren Zustand nicht findet. Ferner werden einige 
diesbeziigliche Experimente besprochen. 


Summary 


Experience, supported by theory, shows that the amount of energy radiated into a room (and, 
therefore, the ways of penetration of outside noise) cannot be deduced from the knowledge of the 
mean amplitude of vibration of the partition. In fact, flexural waves in a partition do not radiate 
unless their velocity be greater than that of sound in air; the “radiation index” of a partition can 
thus be widely different for forced and for free vibrations. A theory of the suspended ceiling is 
attempted, and some observations explained. Finally it is shown that transient phenomena may 
have peculiarities not found in the steady-state, and some relevant experiments are discussed. 


Table de Matiéres dispensable 4 la suite de l’examen et de la discussion 

I. Introduction. de nombreux résultats de mesures de transmission 
II. Notations et définitions. des bruits dans de grands immeubles. 

III. Uné observation concernant l’indice de Dans des batiments modernes ot la transmission 

rayonnement d’un plancher. indirecte est trés importante et dans lesquels le bruit 


se propage & longue distance [1] la question se pose 
de déterminer quelle est, dans une salle ot le bruit 
pénétre, la proportion d’énergie rayonnée par les 
diverses parois. Quand ce ne serait que pour savoir 
ou est le défaut de la construction, quelle paroi doit 
étre modifiée pour améliorer la situation. 

Aucune méthode n’existant actuellement pour la 
mesure des flux acoustiques nous avons entrepris, 
, toutes les fois que nous en avions le temps, de mesu- 
3 I. Introduction rer l’amplitude vibratoire moyenne des parois des 

Nous nous proposons ici d’examiner un certain salles dans lesquelles nous déterminions les ni- 
nombre de problémes dont |’étude nous a paru in- veaux de bruit. 


IV. Calcul du rayonnement d’une paroi plane. 
V. Modes de vibration des parois dans un bati- 
ment. : 
VI. Esquisse d’une théorie du couplage d’une 
paroi légére 4 une paroi lourde. 
VII. Phénoménes transitoires. 
VIII. Conclusions. 
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Nous avons ainsi observé qu’une méme paroi 
peut, a égalité d’amplitude vibratoire moyenne, 
rayonner beaucoup ou trés peu d’énergie selon 
qu’elle fonctionne en transmission directe, indirecte, 
ou vibre sous l’effet d’impacts. 

Avant d’en donner un exemple dont l’exposé 
nécessitera quelques explications préliminaires, si- 
gnalons un autre résultat qui nous a surpris. 

Lors de la reconstruction de l’ilét 4 4 Orléans, 
Varchitecte, M. ABRAHAM, nous avait demandé de 
lui donner les indications nécessaires a la réalisation 
de plafonds suspendus. Ces plafonds furent exécutées 
au moyen de plaques de platre expansé accrochées 4 
des suspentes souples: fil métallique formant une 
boucle pour lui donner lélasticité nécessaire. La 
souplesse des suspentes et leur répartition furent 
choisies de maniére a ce que l’ensemble suspendu ait 
une fréquence propre de l’ordre de 16 Hz. C’était 
done un montage conforme aux données classiques 
en pareil cas. 

Or il se trouva qu’é Vexécution les suspentes 
souples furent, dans certaines piéces, remplacées par 
des fils rectilignes, donc raides. 


Descendants 


! 


4096 


32 64 WA Np NV 
——> Fréquences ¢ 47} 


1024 2048 


Fig. 1. Indice d’isolement de deux planchers pour des sons 
aériens émis dans la salle supérieure ; 
Suspentes souples, 
tere eee Suspentes rigides. 


Ascendants 


Se 64 128 256 512 
Frequences CHz] 


1024 


2048 4096 


Fig. 2. Indice d’isolement de deux planchers pour des sons 
aériens émis dans la salle inférieure; 
Suspentes souples, 
------- Suspentes rigides. 


A notre grand étonnement les mesures mon- 
trérent (Fig. 1 et 2) que, sauf pour les sons trés graves 
(fréquences inférieures 4 128 Hz), les deux struc- 
tures étaient équivalentes. Pourquoi ? 

Un examen théorique des conditions de rayon- 
nement d’un plafond suspendu nous montrera que 
ce résultat n’est pas absurde. 

Ces faits, et d’autres, nous ont done conduit a 
procéder & une étude plus approfondie des problémes 
de transmission des bruits dans le bAtiment en consi- 
dérant les trois étapes suivantes: 


1. Excitation de la structure du batiment par 
un bruit aérien ou un choc. 

2.Transmission interne des vibrations so- 
lides d’un point 4 un autre du batiment. 

3. Dans chaque salle ot le bruit pénétre, rayon- 
nement d’ondes aériennes par les parois qui 
vibrent. 

L’étude expérimentale de ces trois phases du phé- 
nomeéne nous parait indispensable si, au lieu de se 
contenter de constater un état de fait: tel isolement 
a été observé entre telle et telle salle, on veut essayer 
de discerner pour quoi ce résultat est obtenu, et par 
suite comment construire pour obtenir un résultat 
meilleur. 


II. Notations et définitions 


1. Une lettre: X, représente une quantité mesurée 
dans lair. 

La méme lettre affectée d’une apostrophe: X’, 
représente la méme quantité mesurée sur une paroi 
et concerne alors la vibratiou'de cette paroi. 


2.Le niveau sonore aérien est, conformé- 
ment aux normes: 


N = 20log p/Po; (1) 
Po = 2-10~* baryes. 


3.Le niveau de vibration d’une paroi, 
est défini par: 
N’=20log v’/v, (2) 


v’ étant la vitesse vibratoire efficace de la surface 
de la paroi normalement & cette surface, et 

Vp=Po/2 pc=2,4 -10—6 cm/s. (3) 
Le choix de cette valeur de référence a été déter- 
miné par les considérations suivantes : Si nous consi- 
dérons la vitesse vibratoire v d’une particule fluide 
dans un champ acoustique isotrope et diffus, sa 
composante u suivant une direction donnée quel- 
conque est, en moyenne 


u=p/2ee (4) 
p étant la pression acoustique moyenne dans le 
champ acoustique. 
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Si done, dans un tel champ, on mesure la compo- 
sante u (p. e. avec un microphone 4 ruban) le niveau 
sonore sera 

N= 20logu/vp, 


~v» ayant la valeur indiquée plus haut. 


4. Coefficient et indice de rayonnement 
G. Imaginons dans un champ isotrope et diffus un 
plan indéfini P perpendiculaire 4 l’axe Ox (Fig.3); 


Ja moitié des ondes traverse ce plan de gauche a 


droite, autre moitié de droite 4 gauche. Suppri- 
mons ces derniéres ondes. Le champ qui subsiste 
comporte done uniquement des ondes venant de la 
région E (x< 0) et traversant le plan P pour péné- 
trer dans la région R (x > 0). 


Fig. 3. 


Tl est aisé de voir que dans ce champ la relation (4) 
est conservée. Cela résulte des hypothéses de dif- 
fusion et d’isotropie. 

Matérialisons le plan P en lui imposant des vibra- 
tions de flexion telles qu’en chacun de ses points la 
vitesse vibratoire normale 4 P soit égale 4 la com- 
posante selon Ox de la vitesse vibratoire préalable- 
ment existante dans le champ. 

L’espace E peut alors étre considéré comme vir- 
tuel et étre supprimé. 

Dans l’espace R rien n’est changé. Le flux d’éner- 
gie rayonné par l’unité d’aire de P dans l’espace R 
est alors 

=P?/2ec=2ecu”, (5) 


u’ étant la valeur moyenne de la vitesse vibratoire 
de la paroi P. 

Si maintenant nous modifions le mode de vibra- 
tion du plan P tout en maintenant la valeur de uv’, 
le flux rayonné changera et prendra une certaine 
valeur 9,. 


Nous nommons coefficient de rayonnement 


-le rapport 


Pr 
Sins 9 oc Fe (6) 


et indice de rayonnement le nombre 
G=10 log g. (7) 


Le coefficient de rayonnement reste égal a 1 tant 
que la vibration de la paroi P est représentable par 
un ensemble d’ondes provenant d’un espace virtuel 
E dans lequel régne un champ isotrope diffus, 
c.a.d. tant que la distribution spatiale, en ampli- 
tude et en phase, des vitesses vibratoires de P 
correspond & une distribution possible des compo- 
santes selon Ox des vitesses vibratoires dans un 
champ isotrope diffus. 

Dans tous les autres cas g pourra prendre des 
valeurs différentes de 1; G ne sera plus nul. 

Mesure de l’indice de rayonnement en 
salle réverbérante. Lorsqu’une seule paroi d’aire 
S rayonne du bruit dans une salle, ce qui est le cas 
au laboratoire, et aussi entre les salles contigués 
d’un batiment ot la transmission indirecte est né- 
gligeable, la puissance acoustique émise par la pa- 
roi est 

W= Sor. (8) 


D’autre part si le champ acoustique dans la salle 
de mesure est & peu prés isotrope et diffus la pres- 
sion acoustique moyenne p est 


p’=4ocW/A, (9) 
A étant l’absorption totale de la salle. On en tire 
WT p A 


Ta OAT Cia Sin yo) 
9, restant défini par |’équation (5) on trouve 
POS oR EN es 
go ee (525) Sie (Lu) 


N étant alors le niveau sonore aérien mesuré dans 
la salle, N’ le niveau vibratoire moyen de la face 
rayonnante de la paroi 


G =N—N’—3 +10 log A/S. (12) 

5. Indice de transmission interne F, me- 
suré entre deux points 1 et 2 de la structure du 
bAtiment, dont les vitesses vibratoires sont v1’ et v9’ 


F=20logv’,/v’,=N’,—N’,. (13) 


Sa mesure est immédiate. 


6.Indice d’excitation EH. Un champ sonore 
aérien de niveau N crée sur une paroi qu’il touche 
un régime vibratoire de niveau N’. 
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Nous nommerons indice d’excitation le nombre 
E=wN’—N-+3. (14) 
Sa mesure est immédiate. 


7. Décomposition du phénoméne de trans- 
mission directe. 


N, est le niveau sonore aérien dans la salle ou le 
bruit est émis, 

N’, est le niveau moyen de vibration de la face de 
la paroi exposée au bruit dans la salle d’émis- 
sion, 

N’, est le niveau moyen de vibration de la face de 
la paroi qui rayonne du bruit dans la salle de 
réception, 

N, est le niveau sonore aérien dans la salle de ré- 

ception. 

L’indice d’isolement Q de la cloison usuellement 

mesuré est 


Q=N,—N,—10 log A/S. (15) 

On voit de suite que l’on a 
—Q=E+F+6G. (16) 
On a done bien décomposé le phénoméne en ses trois 


étapes. 

La mesure des indices EF, F, G est trés instructive 
et donne des résultats souvent difficiles 4 inter- 
préter. 

Dans le cas des plafonds de l’ilét 4 4 Orléans on 
s’aper¢goit que, pour des isolements équivalents lors- 
que le bruit traverse le plafond de haut en bas ou 
de bas en haut, le mécanisme de |’affaiblissement est 
trés différent. 


Exemple: 


Dans la bande 128—256 Hz on obtient les valeurs 
suivantes 
E F G Q 


de haut en bas sO lL. Date. p 


de basen haut ....... SS Sp) ey ae 
dans la bande 1024—2048 Hz 

E F G Q 
de haut en bas ....... 48 ea Or ea G 
de basen haut ....... mea vtieeh yan | Stes arn is 


Il est frappant de voir les écarts observés dans 
les divers indices, écarts qui montrent que la locali- 
sation des divers affaiblissements est trés différente 
selon le sens de la propagation. 


III. Une observation concernant V’indice de rayon- 
nement d’un plancher. 


Il s’agit d’un plancher formé de poutrelles de 
béton armé coulées sur des corps creux en céramique 
et recouvert d’un parquet en bois sans joints, struc- 


ture trés fréquemment utilisée en France dans les 
batiments & ossature en béton armé. Ce plancher 
pése environ 400 kg/m?. Son indice d’isolement 
moyen entre 100 et 3000 Hz doit done étre de l’ordre 
de 51 db. 

Le batiment ou il est installé (ossature béton, 
remplissages en brique) comporte une transmission 
indirecte importante. 

L’affaiblissemént observé entre une salle du 
2. étage et la salle située au dessous est celui que 
Von aurait dans un batiment sans transmission in- 
directe avec un plancher ayant un indice d’isole- 
ment de 35 db. La transmission indirecte fait done 
perdre 16 db, chiffre énorme. 

D’autre part dans la salle réceptrice le niveau 
vibratoire moyen des six faces: plafond, parois ver- 
ticales, plancher, est a trés peu prés le méme. 

Si Pindice de rayonnement G des six faces était 
le méme, la formule (12) permettrait de le calculer. 
Il apparaitrait 4 peu prés constant entre 100 et 
5000 Hz. Sa valeur serait: + 2 db. 

D’autre part lorsque le plafond est excité par un 
marteau automatique ou un vibromoteur il semble 
y avoir trés peu de transmission indirecte. Nous n’en 
avons pas la preuve expérimentale directe, n’ayant 
pas eu le temps de mesurer la vibration des parois 
verticales. Cependant on observe en général que des 
chocs verticaux appliqués sur le plancher excitent 
peu les vibrations de flexion des parois verticales. 
De plus Vobservation auditive montrait que l’on 
localisait assez bien la position de la source vibra- 
toire située dans la piéce supérieure. Si nous caleu- 
lons alors indice de rayonnement du plafond en 
faisant ’hypothése qu’il est & peu prés seul 4 rayon- 
ner, nous aurons une valeur par excés de cet 
indice. Or le résultat obtenu, en comparant le niveau 
sonore NV mesuré dans la salle au niveau de vibration 
N’ mesuré au plafond, donne les valeurs suivantes: 
de 128 & 512 Hz V’indice de rayonnement G serait 
d’environ — 33 db; au dessus de 512 Hz il croit lé- 
gérement et atteindrait la valeur — 20 dba 4000 Hz. 

Que conclure de ces deux observations ? 

Pour ]’émission aérienne, si les six faces avaient 
le méme indice de rayonnement l’isolement normal 
du plancher n’aurait été réduit que de 8 db. Or la 
réduction peut étre évaluée 4 environ 16 db. 

Il faut donc admettre l’une des deux hypothéses 
suivantes : 
ou bien le plancher est exceptionnelle- 
ment mauvais eu égard a son poids et rayonne 
& peu prés autant que les autres parois (¢ = + 2). 
Mais alors, sous l’effet d’impacts ou de 
vibrations son indice de rayonnement est 
de 30a 25 db plus faible, selon la fréquence. 


Pa eo 
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Ou bien le plancher a un indice d’isole- 
ment normal,dten grande partie a ce qu’ilrayon- 
ne peu, puisque dans ce cas son indice de 
rayonnement doit étre trés inférieur a ce- 
lui des cloisons qui, sans avoir de plus grandes 
amplitudes vibratoires, sont responsables de la 
perte d’isolement global de 16 db. D’ot vient alors 
cette différence entre le rayonnement des cloisons 
et celui du plancher ? 
| Dans l’une comme dans lautre hypothése la que- 
stion se pose: comment rayonne une paroi? Peut 
elle, selon les cas, avoir des indices de rayon- 
nements trés différents ? 


IV. Caleul du rayonnement d’une paroi plane. 


Considérons un plan P indéfini (x =0) parcouru 
par des ondes de flexion qui se propagent vers les y 
croissants & la vitesse ¢,. 

La vitesse vibratoire des points de ce plan est di- 
rigée suivant Ox et a pour valeur 


u’=u’,-exp [iw(¢—y/¢)]. (1) 
Dans le demi espace x > 0, la vibration du plan P 
crée un champ acoustique dont le potentiel des vi- 
tesses ® peut s’écrire 
® =e Y (x,y). (2) 
L’équation générale du son s’écrit alors 
ey egy a? 
ox? ay? 4 at alls 


c étant la vitesse de propagation du son dans lair. 
Sur le plan x =0, la vitesse vibratoire est imposée 
par la condition (1). On doit done avoir pour « = 0 
oY 
0x 


(3) 


=u',-exp(—ioy/4). (4) 


Ecrivons V’ sous la forme 

Y= A-exp[i(ax +8 y)]. (5) 
En portant cette valeur dans les relations (3) et (4), 
il vient 

a +p? = alle? (6) 
et 
in A-exp(iPy) =u’, -exp(—ioy/c,). 

Cette derniére relation doit étre satisfaite quel que 
soit y. On en tire donc 


8 = — a/c, (7) 
valeur qui, portée dans (6) donne 
1 1 
nese ot 
a2 = © ( 2 2 (8) 
On a enfin 
= —iu’,/e. (9) 


Deux cas peuvent se présenter: 

1. c, > c. La vitesse de propagation des ondes de 
flexion est supérieure a la vitesse du son dans Lair. 
« est alors réel et peut s’écrire 


a) 


Kania cos 0 (10) 
avec sin 0 =c/e¢,. 
On a alors 
: & 
rn iseos en (11) 
VY=A.exp|—i : (x cos0+ysin0)}. (12) 


Le champ acoustique est constitué par une onde 
plane progressive dont la direction de propagation 
fait avec Ox l’angle 9. 

Le flux rayonné par la paroi est alors: 


ue pcuy = } 
?r = 2cos0 2cosd Me) 


ou o, =2ecu’,?/2, u’,/V 2 étant la valeur efficace de 
1 1 1 

la vitesse wu’. 
Le coefficient de rayonnement est 


g =1/2cos 0 (14) 
Pour) cy) = ce 6=0 g=4, 
quand c, diminue, 0 croit, g augmente 
Pour, ¢,=2e/) 3 §=n7/3 g=1 
lorsque c, —>¢ 0 > 7/2 g-> 


2.c,<c. « est alors un imaginaire pur et peut 
s’écrire 


sinh (15) 


. @ 
== 
c 


avec 
cosh =c/e,. 
Il vient alors 
Cc 


, 
A =— -UW4, 


Eeeinhn (16) 


Y= A .exp —" (x sinh ) + iy cosh ¥) |. (17) 


L’onde émise par le plan P est plane et non uni- 
forme (Fig. 4). La surface d’onde est paralléle au 
plan xOz. L’onde progresse dans le sens des y crois- 
sants & la vitesse c,. Les particules aériennes ont des 
trajectoires elliptiques dont le grand axe est paral- 
léle & Oy et le petit axe paralléle 4 Ox. Le sens de 
parcours est indiqué sur la figure. La composante 
de la vitesse dans le sens Oy est en phase avec la 
pression. La composante dans le sens Ox est en 
quadrature avec la pression. Le rayonnement 
de la paroi est done nul. 
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Ces ondes planes et non uniformes, dont l’ampli- 
tude décroit exponentiellement vers les x croissants, 
se rencontrent dans la théorie de la réflexion totale. 
En acoustique on les trouve également lorsqu’on 
étudie la propagation du son entre deux parois ab- 
sorbantes [2] ou la décomposition en 8 ondes planes 
obliques d’un mode de vibration propre amorti d’un 
parallélépipéde a parois absorbantes. 
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Ul \ 
\ 
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\_ 7 Surface d’ onde - 
maxima ‘ 


~ 


Sens de propa - 
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Fig. 4. Onde engendrée par une paroi que parcourt une onde 
de flexion dont la vitesse est inférieure & la vitesse 
du son dans lair. 


Il s’agit 14 d’un phénoméne qui ressemble a la 
houle, c.a.d. qui reste en surface et pénétre peu 
dans la profondeur du fluide. I différe cependant 
de la houle trochoidale car cette derniére n’est pas 
un mouvement irrotationnel comme l’est notre 
onde. 


Affaiblissement de l’onde dans le sens Ox. 


Il convient d’examiner si l’affaiblissement du 
champ sonore dans le sens des x croissants est assez 
rapide pour que, dans la pratique, ce phénoméne 
ait de importance. Cet affaiblissement est donné 


par le facteur exp (— “= sinh ’). Il est done, en 
décibels, 
D= 8,7" sinh ¢ = 54,7 © sinh. 


Pour avoir un affaiblissement de 10 db il faut done 
s’éloigner 4 une distance de la paroi égale & 
OBS: 
“10 ~ 5,47 sinh ¢ * 
Le tableau suivant donne quelques valeurs de x19/4 
en fonction de ¢,/c. 


Chi clestie 1 
BAIS oak 9 


0,99 0,95 0,9 0,5 
1,29 0,568 0,38 0,106 
On voit que, méme pour des vitesses de propaga- 
tion c, trés peu inférieures & la vitesse du son dans 
Vair, l’affaiblissement est extrémement rapide. 


Cas d’ondes stationnaires. 


Dans la pratique, nous n’aurons pas affaire 4 une 
onde de flexion progressive le long d’un plan in- 
défini, mais 4 un systéme d’ondes stationnaires dans 
une paroi limitée. Ces ondes stationnaires sont re- 
présentables par la superposition d’ondes progres- 
sives. Si ces ondes progressives ont une vitesse de 
propagation c, supérieure a c, la paroi rayonne; si ¢, 
est inférieur 4 c, la paroi ne rayonne pas. 

En fait la paroi se trouve découpée en mailles par 
des lignes nodales (Fig. 5). D’une maille aux mailles 
contigués la vibration de la paroi est en opposition 
de phase: lorsqu’une maille est en saillie, ses voi- 
sines sont en creux. 

Lorsque la vitesse de propagation ¢, des ondes 
de flexion dans la paroi devient inférieure 4 c¢, les 
mailles sont petites. L’air chassé par une bosse, au 


‘A 
A ( 
Plan Coupe AA 


Fig. 5. Mouvements de l’air au voisinage d’une paroi sou- 
mise & des ondes stationnaires de flexion dont les 
lignes nodales sont trop rapprochées pour que la 
paroi rayonne. 


| 


| 
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lieu d’engendrer une onde progressive, s’écoule 
simplement dans les creux voisins. C’est exactement 
le méme phénoméne que |’on observe dans les fré- 
quences graves avec un haut-parleur sans écran. 
Il ne rayonne plus: l’air chassé devant la membrane 
la contourne et s’écoule derriére elle. 


Fig. 6. Si les haut-parleurs blancs et noirs sont en phase, 
le systéme rayonne quelle que soit la fréquence. 
S’ils sont en opposition de phase, le rayonnement 
n’apparait que pour les sons aigiis, dont la Jongueur 
d’onde dans Vair est inférieure a 2b. 


Une expérience simple permet d’entendre ce phé- 
noméne. La paroi est schématisée par un écran sur 
lequel on monte des petits haut-parleurs disposés 
en damier (Fig. 6). On alimente séparément les deux 
groupes figurés en noir et blanc. Un commutateur 
inverseur permet de connecter ces groupes soit en 
phase, soit en opposition. Dans le premier cas la 
paroi fonctionne dans l’ensemble comme un piston 
et rayonne de l’énergie. Dans le second cas le 
systéme figure une paroi dans laquelle existe un 
systéme d’ondes stationnaires dont les mailles sont 
carrées, les lignes nodales passant entre les haut- 
parleurs. Le cété 6 des carrés ainsi définis est alors 
égal a ,/2; A, étant la longueur d’onde des ondes de 
flexion dans la paroi. Les mailles du systéme ont 
alors les dimensions qu’elles auraient si la vitesse 
de propagation des ondes de flexion était c, = A,n 
=2bn, n étant la fréquence. La coupure se pro- 
duira done & la fréquence ny =c/2b pour laquelle 
la longueur d’onde 2 du son dans l’air est A) = 2b. 

Quand la fréquence n est inférieure a mo, le son 
produit est beaucoup plus faible lorsque les groupes 
de haut-parleurs sont en opposition que lorsqu’ils 
sont en phase. Quand n est supérieur a no, l’inten- 


sité rayonnée dans les deux cas est sensiblement la 
méme. 

Ce montage peut d’ailleurs étre utilisé comme 
émetteur passe-haut, la fréquence de coupure étant 
réglée par la distance 6 entre haut-parleurs. 

Une derniére remarque doit encore étre faite: 
dans la pratique, la présence d’obstacles plus ou 
moins absorbants au voisinage de la paroi déforme 
le champ acoustique théorique. Les ondes aérien- 
nes produites deviennent légérement obliques, 
faisant ainsi apparaitre le rayonnement nécessaire 
pour remplacer l’énergie absorbée. 


V. Modes de vibration des parois dans un batiment 


Pour voir comment le rayonnement d’une paroi 
réelle peut se modifier selon les circonstances il faut 
examiner comment les parois d’un batiment sont 
susceptibles de vibrer. 


Vibrations forcées. 


Ce cas est celui d’une paroi soumise sur toute sa 
surface & des pressions acoustiques, comme cela se 
produit dans la transmission directe d’ondes 
aériennes par une paroi simple. La distribution spa- 
tiale des amplitudes et des phases de la vibration est 
alors imposée par celle qui existe dans l’air sur la 
face exposée au bruit. 

La vitesse de propagation des ondes de flexion 
libres dans la paroi n’intervient alors que pour régir 
dans une certaine mesure l’amplitude des vibrations 
de la paroi, ainsi que l’a montré CREMER, [3]. 

A chaque onde incidente correspondra une onde 
transmise ayant la méme direction de propagation 
mais dont l’amplitude sera fonction de l’angle d’in- 
cidence 8. Le champ acoustique rayonné par la 
paroi restera dans une certaine mesure diffus, mais 
ne sera plus isotrope. Dans ce champ transmis 
le flux rayonné par la paroi n’est done plus lié a la 
vitesse vibratoire moyenne de la paroi par la rela- 
tion 9 =2pcu”, qui postule lisotropie. 

Le coefficient de rayonnement g pourra done étre 
différent de 1, l’indice G différent de 0. 


Vibrations libres 


Dans le cas de vibrations libres les mouvements 
de la paroi résultent de la superposition d’ondes de 
flexion progressives qui cheminent a la vitesse c¢,. 

Ce sera le cas d’une paroi simple attaquée par un 
marteau ou une source vibratoire localisée. 

La vitesse de propagation c, des ondes de flexion 
est, on le sait, proportionnelle a la racine carrée de 
la fréquence. On peut donc l’écrire 


ee 
c, =cyn/ng, 


ny étant la fréquence pour laquelle c, = c. 
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Par suite, pour les vibrations dont la fréquence 
est supérieure 4 mo, la paroi rayonnera. Pour celles 
dont la fréquence est inférieure 4 no, elle ne ray- 
onnera pas. 

La fréquence de coupure ny se mesure aisément 
pour peu que l’on ait de Voreille en repérant la hau- 
teur de la fondamentale et des premiers partiels 
dune verge taillée dans un élément de la paroi. 
Nous en donnons quelques valeurs ci-dessous: 
Plaques de platre expansé (plafonds 

de Vilot 4.4 Orléans) Ny = 680 Hz 
Placoplatre de 1 cm d’épaisseur.. nm) =3080 Hz 
Tole d’acier de 1 mm d’épaisseur 

(valeurjcalculée) shies cian 1, = 6400 Hz. 


er 


Bruit produit par une paroi en vibration 


libre. 


Pour nous rendre compte de l’importance du phé- 
noméne il est bon de donner quelques chiffres. Nous 
avons done calculé le niveau sonore N (mesuré au 
moyen d’un microphone-pression) qui existe devant 
une paroi indéfinie parcourue par des ondes de 
flexion dont le niveau vibratoire N’ reste constant. 
Nous donnons dans les tableaux suivants la valeur 
du terme N—W’ en fonction de la fréquence n pour 
une paroi dont la fréquence de coupure ny est égale 
& 1000 Hz, valeur qui est dans le domaine de la pra- 
tique comme le montrent les chiffres donnés plus 
haut. 
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Fig. 7. Indice de rayonnement d’un plancher muni d’un 
plafond suspendu par attaches souples dans le cas 
d’une émission de sons aériens dans la salle supéri- 
eure. 

Lorsque n est supérieur & my ,le niveau N est indé- 
pendant de la distance du microphone 4 la paroi 
(onde progressive). 

Lorsque n est inférieur & ny, N décroit lorsqu’on 
s’éloigne de la paroi. Nous avons caleulé N—N’ 
pour des distances x de 0,50—1—2—4 m, 


n= no n [Hz] N—N’{db] 
oo —6 
10 000 — 5,5 

5 000 —5 

2 000 —3 

ees: 0 

1111 4 

1010 14 

1 000 oo 

nN =n,) n [Hz N—N’ [db] 

&0,50m admis sans a4m 
1000 +o + co + co + co 
990 +6 —2 —18 —850 
909 —19 — 42 —8s8 —180 
800 —32 — 64 —129 —258 
667 —4l —79 —155 —306 
500 —46 — 86 —167 —328 
333 —47 — 85 —160 —38312 
250 —46 —80 —I150 —289 
100 —40 — 64 —112 —208 
50 —36 . —54 — 89 —159 


Ces chiffres, 4 notre avis, n’ont pas besoin de 
commentaires. 


Parois anisotropes 


Dans le cas de parois anisotropes (bois, struc- 
tures comportant des poutrelles, etc. . . .) la vitesse 
de propagation des ondes de flexion dépend de la 
direction de propagation des ondes. On aura done 
deux directions dans lesquelles n, sera maximum ou 
minimum. Entre ces deux valeurs de n, les ondes 
rayonneront ou ne rayonneront pas suivant leur 
direction. La coupure sera donc, en pratique, moins © 
marquée qu’avec une cloison isotrope. 


Vibrations semi-libres 


Dans beaucoup des cas pratiques et principalement 
dans les cas ot: on observe des transmissions indi- 
rectes importantes, les parois susceptibles de 
rayonner sont relativement légéres par rapport 4 
Vossature du batiment. On peut estimer que leur 
présence ne modifie pas fondamentalement les vibra- 
tions de cette ossature. Par suite les lignes ou les 
points d’attache de ces parois 4 l’ossature sont en 
vibration forcée, ou, tout au moins, subissent un 
régime vibratoire dont la répartition spatiale dépend 
beaucoup plus des propriétés de l’ossature que de 
celles de la paroi. 

La paroi prendra done un régime de vibration que 
nous hommerons «semi-libre», dans lequel ses pro- — 
priétés joueront un réle d’autant plus important 
que les liaisons avec |’ossature seront moins nom- 
breuses et moins rigides. 


A ce point de vue une liaison par points isolés 
sera moindre qu’une liaison par lignes; un appui 
moindre qu’un encastrement. 

Dans le cas des plafonds suspendus de lilét 4, 

A Orléans, dont nous avons parlé au début de cet 
article — suspentes rigides dans un cas, souples 
dans l’autre, mais réparties de la méme maniére — 
nous avons tracé les diagrammes de l’indice de rayon- 
nement des plafonds. Ces diagrammes sont va- 
lables, toutes les mesures faites dans ce batiment 
ayant montré que la transmission indirecte y est 
faible. Ils sont, dans l’ensemble, assez confus et 
difficiles & interpréter mais tous montrent une 

- augmentation brusque de l’indice de rayonnement 
(15 & 30 db) lorsque l’on passe des octaves inférieurs 
a 1024 Hz aux octaves supérieurs. Or, nous l’avons 
indiqué plus haut, la fréquence de coupure np des 
plaques de platre employées pour ces plafonds est 
de 680 Hz. 

Nous donnons deux de ces diagrammes correspon- 
dant l'un (Fig. 7) 4 un plancher & suspentes souples 
dans le cas d’une émission aérienne dans la salle 
supérieure, l’autre (Fig. 8) 4 un plancher a suspentes 
rigides dans le cas d’impacts par un marteau. . 

Il semble donc bien que le fait de l’accrochage du 
plafond par points isolés ait été plus important que 
la nature des suspentes. Ce mode d’accrochage per- 
met en effet au plafond de prendre des déformations 
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% Fig. 8. Indice de rayonnement d’un plancher muni d’un 
t plafond suspendu par attaches rigides dans le cas 
dimpacts. ? 
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aN Fig. 9. Schéma d’un plancher avec plafond suspendu. 


J. BRILLOUIN: PROBLEMES DE RAYONNEMENT EN ACOUSTIQUE DU BATIMENT 73 


presque libres méme en cas d’excitation par bruit 


-aérien, c.a.d. de vibration forcée du plancher: 


Vexistence de la coupure dans ce cas en est la preuve. 


VI. Esquisse d’une théorie du couplage d’une paroi 
légére a une paroi lourde 

Ce probléme est l’extension de celui de l’accro- 
chage d’un plafond. Nous nous proposons d’exami- 
ner comment une liaison par lignes ou par points 
isolés transmet les vibrations d’une face a l’autre de 
la paroi composite. 

Pour simplifier nous considérerons d’abord un 
probléme 4 deux dimensions. Soit A (Fig. 9) un 
plancher porteur relativement massif auquel un 
plafond B relativement léger est accroché par 
une série de suspentes L réguliérement espacées. 
Nous supposons ces suspentes transversalement as- 
sez souples pour ne transmettre au plafond que des 
mouvements verticaux. 

La différence de masse entre les deux éléments 
nous permet d’admettre, en premiére approxima- 
tion, que les vibrations du plancher A ne sont pas 
sensiblement modifiées par la présence du plafond. 

Lorsque Je plancher A atteint par des vibrations 
est parcouru par un onde de flexion progressive se 
propageant vers les x croissants, l’élongation ver- 
ticale z d’un point de A est de la forme 


2=A -exp [iw(t—x/e,)]. (1) 


La distance entre les points d’attache L étant a, 
leurs abscisses sont: x,,—= ma (m entier, positif, 
négatif ou nul). 

Leurs mouvements verticaux 2, sont alors 


2m =A-exp [ia(t—ma/c,)]. (2) 


Effet de la liaison 


Utilisons d’abord un mode de raisonnement au- 
jourd’hui classique dans les problémes de propaga- 
tion d’ondes dans les structures périodiques. Il con- 
siste 4 remarquer que la vitesse de propagation c, 
de l’onde du plancher n’est que partiellement con- 
nue par les effets de cette onde sur le plafond. 

En effet, le déplacement z,, d’une suspente ne 
change pas, si l’on remplace la vitesse c, par une 
vitesse c’, liée A la précédente par la relation 


aw aw 


exis + 2kr, 
Cy Cy 
c’est-a-dire 
1 1 k 
aS + (3) 
Cy (ia an 
(n = w/27). 


Ceci revient 4 dire que l’on peut, sans que le pla- 
fond s’en apergoive, remplacer une onde ayant la 
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vitesse c, et la longueur d’onde A, par une autre, 
de vitesse c’, et de longueur d’onde 2‘;, avec 


ee (4) 

C’est ce que fait comprendre la Figure 10. Il y 

apparait avec évidence que l’onde longue (trait in- 

terrompu) se déplacgant rapidement produira le 

méme effet sur les points L que londe courte (trait 

continu) de méme fréquence qui se propage 9 fois 
plus lentement. 


Fig. 10. Les suspentes L du plafond prennent le méme 
mouvement lorsque le plancher est parcouru par 
les deux ondes | (trait continu) cheminant lente- 
ment, et 2 (trait interrompu) cheminant neuf fois 
plus vite; 

Vitesse 1 cm/s, 

were eee Vitesse 9 cm/s. 


Par suite, le plafond est inapte 4 déceler le sens 
de propagation et la valeur exacte de la vitesse de 
propagation de l’onde du plancher. Nous pouvons 
done simplifier notre étude en prenant pour c’, la 
vitesse la plus grande en valeur absolue qui soit 
compatible avec la relation (3), c.a.d. celle qui 
correspond a la longueur d’onde 


Kj 2a. 


Par conséquent tout se passe dans le 
plafond comme si la vitesse de propaga- 
tion des ‘ondes dans le plancher était tou- 
jours supérieure a 2an. 

Cette vitesse apparente ne saurait done étre in- 
férieure & la vitesse du son dans l’air que pour les 
fréquences n <c/2a, c’est-d-dire pour les ondes dont 
la longueur A dans lair satisfait & 4> 2a. 


Forme des ondes dans le plafond 


Comme, dans cette esquisse de théorie, nous con- 
sidérons un phénoméne illimité (de x=—o Aa 
x==-+ co) dans une structure périodique réguliére, 
les ondes du plafond doivent étre représentables 
par la superposition d’ondes sinusoidales se pro- 
pageant sans affaiblissement. Seules les ondes 
ayant des vitesses de propagation c’, régies par la 
relation (3) 


satisfont aux conditions de périodicité temporelle et 


spatiale imposées par le mouvement des suspentes, 
savoir: 
2m = Ay - exp [iw(t— ma/c’;)]. 

Or les fonctions f; = exp (iw x/c”,) forment, dans 
Vintervalle (0, a) qui sépare deux suspentes, une 
suite de fonctions orthogonales, car on a: 

a 


[exp (—iase’,)exp (iox/e'n)dix=0 
J 


lorsque c’, et c’,, satisfont a la relation (8). 
Pour justifier en toute rigueur la représentation 
du mouvement du plafond par une série 


+ 00 
exp (iw?) yr -exp (—iwx/e’;) 


k =— 0 


il faudrait démontrer que la suite des fonctions f;, est 
compléte et que la série est convergente dans l’inter- 
valle considéré. Faute de cette démonstration nous 
nous contenterons dans la présente esquisse de la 
présomption suivante: lorsque les mouvements des 
suspentes sont en phase, la vitesse ¢, est infinie, les 
fonctions fx prennent la forme f, = exp (i2knx/a). 
Notre série se réduit done & une série de Fourier qui, 
elle, jouit des propriétés nécessaires. Nous ad- 
mettrons donc que notre représentation reste va- 
lable dans les autres cas. 

Les vitesses de propagation des diverses ondes 
ainsi séparées sont en nombre infini. La plus grande 
est, en valeur absolue, supérieure a an/2. 

Il n’y aura done qu’un nombre limité 
de ces ondes qui aient une vitesse de pro- 
pagation supérieure 4 la vitesse du son 
dans l’air, et par suite qui soient susceptibles de 
fournir un rayonnement dans l’air situé sous le pla- 
fond. Au contraire, toutes les ondes (théoriquement 
en nombre infini) dont les vitesses de propagation 
sont inférieures ala vitesse du son dans |’air, nerayon- 
neront pas d’énergie ou en rayonneront trés peu. 

Il est d’autre part intuitivement évident que les 
ondes qui auront la plus grande amplitude seront 
celles dont la vitesse de propagation ec’; sera le plus 
proche de la vitesse de propagation des ondes de 
flexion libres du plafond. Ce phénoméne & la fois 
spatial et temporel est en effet analogue a une ré- 
sonnance. 

Par suite ilaura interét a ce que la paroi 
constituant le plafond ait la fréquence de 
coupure n la plus élevée possible. 

D’autre part, la série des vitesses compatibles 
avec la relation (3) sera d’autant plus serrée que 
Vécart a entre suspentes sera plus grand. Done, plus 
cet écart sera grand, plus on aura de chance qu’une 
onde ayant une vitesse trés proche de la vitesse des 
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ondes propres du plafond soit excitée de préférence 
& des ondes plus rapides. Cette considération rejoint 
une évidence intuitive: plus les suspentes seront 
écartées moins elles seront nombreuses; moins il y 


-a de contacts entre le plafond et le plancher plus 


la vibration du plafond se rapproche de la forme 
d’une vibration libre qui, si la fréquence est in- 
férieure 4 my), ne rayonne pas. 

Il y a done intérét & espacer les sus- 
pentes le plus possible. 

Réaliser ces suspentes par points isolés est évi- 
demment nécessaire si l’on veut que les avantages 
mis en lumiére ci-dessus existent pour des ondes 
cheminant dans n’importe quelle direction dans le 
plancher. 

L’analyse sommaire des phénoménes que nous 
venons de donner permet donc de comprendre les 


_faits observées dans les planchers de l’Il6t 4 a 


Orléans. 

Elle nous indique également une voie: Toutes les 
fois que les conditions statiques d’une construction 
imposent & une paroi une certaine raideur: murs 
portants par exemple, cette paroi aura une fré- 
quence de coupure relativement basse (on trouve 
dans l’article de Cremur déja cité des données pré- 
cieuses 4 ce sujet). Il sera donc avantageux de mas- 
quer cette paroi par une paroi plus mince, attachée 
par points isolés 4 la paroi portante, et qui, elle, 
sera susceptible de vibrer sans faire de bruit. 

La technique actuelle des murs composites qui 
jusqu’ici ont surtout été étudiés au point de vue de 
V’étanchéité et de l’isolation thermique, parait se 
préter aisément a de telles réalisations. Des mesures 
en cours dans un batiment construit de la sorte nous 
ont déji montré que de tels murs — bien que n’ay- 
ant pas été spécialement étudiés au point de vue 
acoustique — paraissent doués de qualités intéres- 
santes. 


VII. Phénoménes transitoires 


Nous avons jusqu’ici examiné le cas de phéno- 
ménes permanents et nous avons montré que, au 
voisinage d’une paroi qui vibrait sans rayonner, il 
existe un champ acoustique dont l’intensité décroit 
exponentiellement lorsqu’on s’éloigne de la paroi. Il 
s’agit done la d’un champ acoustique localisé au 
voisinage de la source, phénoméne général que |’on 
rencontre dans |’étude de toutes les sources sonores. 
Un tel champ est stationnaire. Son entretien de- 
mande uniquement que le mouvement de la source 
se perpétue, mais n’emprunte aucune énergie 4 la 
source. L’énergie acoustique existant dans ce champ 
est stationnaire, ,,stagnante“. Elle a été mise en 
place en prélevant de ]’énergie sur la source au mo- 


ment ot le mouvement de cette derniére a débuté. 


Elle se dissipera sous forme d’ondes progressives 


lorsque la source s’arrétera. 

Nous retrouvons 14 le probléme que j’ai récem- 
ment traité dans le cas de sources sphériques [4]. 

Tl parait alors certain que les phénoménes que 
lon observera, et leurs effets sur l’audition, pourront 
étre trés différents selon que l’on aura affaire 4 une 
vibration prolongée de la paroi ou 4 des trains d’onde 
brefs comme ceux qui résultent d’impacts. 

Nous n’avons pu faire encore des expériences systé- 
matiques sur cette question, mais pouvons signaler 
le résultat de mesures qui en prouvent |’existence 
et importance pratique. 

Dans l’article visé ci-dessus j’avais indiqué que les 
ondes transitoires de diffraction qui se produisent 
sur les nombreux reliefs des parois de studios du 
type dit «polycylindrique» paraissaient de nature 
& avoir une influence sur la réverbération. Or 
au début de l’été dernier, M. Conturiz, Ingénieur 
au Service des Batiments de la Radiodiffusion 
m’avait signalé certaines singularités observées dans 
les studios du Centre Rop1n en cours d’achévement. 
Ayant alors tenté d’en donner une explication 
théorique — tentative d’ailleurs incompléte et in- 
suffisante — je la soumis & cet Ingénieur qui décida 
de faire exécuter quelques mesures en vue de pré- 
ciser les données expérimentales du probléme. 

Ce travail, effectué avec le plus grand soin par 
M. Borsarp, fut le suivant: Deux studios furent 
choisis, l'un ayant des parois planes, l’autre pourvu 
de nombreux reliefs en forme de calottes sphériques. 
Dans chacun d’eux on enregistra la courbe de ré- 
verbération d’une part aprés émission prolongée 
d’un son de bande filtré, d’autre part en utilisant 
un claquoir. Les enregistrements furent faits par 
bandes d’octaves, de 100 & 12800 Hz. 

Dans le premier studio (parois planes) les mesures 
faites au claquoir donnent dans |’ensemble des du- 
rées de réverbération trés légérement plus longues 
que celles faites en émission soutenue. L’écart est 
faible, et reste de l’ordre de grandeur des écarts 
expérimentaux normaux. Dans le second au con- 
traire (studio bosselé) les durées de réverbération T, 
mesurées au claquoir sont 4 toutes les fréquences 
nettement plus faibles que les durées T, mesurées 
en émission soutenue. Le rapport T,/T, est en moy- 
enne égal 4 0,46; les valeurs extrémes observées sont 
0,30 (bande 100/200 Hz) et 0,61 (bande 1600/3200 
Hz). 

Tl y a donc la un effet systématique et important 
qui peut étre une des causes de la qualité (ou des 
défauts, c’est aux artistes 4 en juger) des studios 
polyeylindriques. 
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Cet effet mériterait d’étre étudié avec soin tant 
théoriquement qu’expérimentalement. 

En tous cas, on doit en tenir compte immédiate- 
ment. Si, dans une salle a parois bosselées les équa- 
tions usuelles de la réverbération ne s’appliquent 
plus d’une maniére générale au phénoméne de la 
décroissance de l’intensité sonore, elles seront égale- 
ment en défaut pour l’évaluation du niveau sonore 
existant en régime permanent, et, en particulier 
pour la comparaison du niveau en régime permanent 
stable 4 celui qui est produit par des impacts répétés. 
Par suite, tant que ce phénoméne n’aura 
pas été éclairci et bien étudié, on doit 
prohiber l’emploi comme salle de mesure 
d’un local contenant un grand nombre 
déléments diffractants, qu’ils soient placés 
' sur les parois ou 4 l’intérieur de la salle. 

De plus, ce coup de sonde dans un domaine encore 
insuffisamment exploré montre avec évidence qu’il 
faut s’attendre 4 observer des phénoménes singu- 
liers, 4 avoir des surprises, lorsque l’on comparera le 
comportement de parois en régime permanent ou 
sous effet de chocs. 


VIII. Conclusions 


Sans pouvoir aboutir a des conclusions défini- 
tives, notre étude conduit cependant a un certain 
nombre de résultats: conseils de prudence, possibi- 
lité de réalisations techniques intéressantes, champ 
de recherches nouvelles, tant dans le domaine du 
laboratoire que de la construction. 


A. Problémes de mesure 


1. wutiliser qu’avec beaucoup de circonspection les 
salles douées d’une diffraction importante. 

2. se méfier des mesures de transmission effectuées 
avec un microphone placé trop prés de la paroi: 
ce procédé autrefois utilisé en acoustique du bati- 
ment est encore parfois employé pour l’étude de 
structures légéres destinées 4 la construction de 
earlingues d’avions. Or les téles ont en général 
une fréquence de coupure trés aigtie. On est done 
fondé & se demander si cette méthode, surtout 
lorsqu’on l’applique & des structures complexes, 
ne risque pas de donner des résultats trompeurs, 
le mesure ayant porté sur le champ stationnaire 


au voisinage immédiat de la paroi et non sur le 
rayonnement de la paroi. 


B. Probleme d’observation 


Il serait fort intéressant que la mesure des vibra- 
tions des parois fut systématiquement effectuée 
dans tous les immeubles ot les techniciens ont l’oe- 
casion de faire des mesures de transmission sonore. 
Ce travail d’«observation clinique» serait & notre 
sens un des plus fructueux: lorsqu’on interroge la 
nature il convient d’étre exigeant et de la pousser 
dans ses derniers retranchements. 


C. Problémes de construction 


1. Méme lorsqu’il s’agit de la transmission directe 
par une paroi simple, et a priori dans tous les 
autres cas, le rayonnement d’une paroi dépend 
essentiellement de la vitesse de propagation des 
ondes libres de flexion dans cette paroi. 

2. La valeur des indices d’isolement mesurés au 
laboratoire en transmission directe ne donne de 
renseignement valable ni sur le comportement 
des parois en transmission indirecte ni sur leur 
rayonnement sous leffet d’impacts. 

3. L’emploi de parois relativement légéres et soup- 
les, ayant une fréquence de coupure mp assez 
élévée, doit permettre de réaliser des structures 
intéressantes (planchers, parois composites). Ces 
parois pourraient étre efficaces sans étre «flot- 
tantes» si elles ne sont reliées 4 Vossature ot a 
élément porteur que par des points isolés et 
assez espacés. ‘ 

(Regu le 28. Mai 1951.) 
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BEMERKUNGEN ZUR GEOMETRISCHEN RAUMAKUSTIK 


Von E. MEYER und W. KUHL 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Die vorliegende Arbeit befaBt sich mit dem Problem der ersten Riickwtrfe im Anschlu8 an 
den direkten Schallstrahl Schallquelle—Zuhérer. Man kann diese Riickwiirfe durch Anbringung 
geeigneter Schallspiegel insbesondere auf den oberen Teilen der Seitenwande erheblich vermehren. 
Betrachtungen tiber die dadurch mégliche Verstarkung der innerhalb der ersten 50 ms ankom- 
menden Schallstrahlen und tiber die Energiedichte im diffusen Schallfeld unter diesen Bedin- 
gungen werden angestellt. 

Im Anschlu8 an die physikalischen Ausfihrungen werden eine Reihe von praktischen Bei- 
spielen besprochen, in denen der Einsatz der Schallspiegel zu bemerkenswerten Verbesserungen 
der Verstandlichkeit von Sprache gefiihrt hat. Es sind dies die wiederaufgebaute Staatsoper 
Hamburg, der grofe Saal im Studentenhaus in Bonn, sowie das Landestheater in Hannover. 


Summary 


The paper deals with the first sound reflections following the direct beam from source to 
listener. The number of these reflections can be greatly increased by sound mirrors suitably tilted 
and placed at the upper parts of the side walls in the neighbourhood of the sound source. If these 
first reflections occur in the 50 ms interval after the arrival of the direct sound a reinforcement 
takes place. Some formulae are given relating this effect to the energy density in the diffuse 
sound field. 

The second part of the paper describes some practical applications to rooms, in which sound 
mirrors improve the hearing conditions and especially the intelligibility of speech. They are: the 
rebuilt opera house at Hamburg, the big hall in the students hostel at Bonn, and the theatre at 
Hannover. 


Sommaire 


La présente note traite le probléme des premiéres réflexions du son qui ont lieu aprés la 
transmission directe source sonore-observateur. On peut augmenter considérablement ces réflexions 
en disposant des miroirs acoustiques appropriés sur les parties supérieures des parois. On étudie 
le renforcement qui peut se produire pendant les 50 premiéres millisecondes et la densité d’énergie 
existant dans le champ sonore diffus dans ces conditions. 

La seconde partie donne une série de réalisations pratiques montrant l’amélioration notable 
de Vintelligibilité que l’on peut obtenir en installant des miroirs acoustiques. Il s’agit de Popéra 
de Hambourg, de la grande salle de la Maison des Etudiants & Bonn ainsi que du théatre de 
Hanovre. 


rift man fertige Raume oder Planungen von 

Raumen auf ihre Horsamkeit hin, so tut man 
dies im allgemeinen nach zwei Gesichtspunkten, und 
zwar hinsichtlich der geometrischen Raumakustik 
und hinsichtlich der statistischen Raumakustik, d.h. 
man untersucht den Verlauf aller Schallstrahlen von 
der Schallquelle zum Beobachter hin, und man be- 
faBt sich mit den Fragen des Nachhalls. Fiir den 
ersten Punkt ist die Raumform das Entscheidende, 
wahrend im zweiten Fall die Schallabsorption der 
Wandbekleidung im Vordergrunde steht. Es gab in 
der Entwicklung der Raumakustik Zeiten, in denen 
man bald die eine, bald die andere Betrachtungs- 
weise starker herausstellte, ja sie zum Teil sogar als 


die allein zweckmaBige ansah und dementsprechend 
bei den praktischen Anwendungen Riickschlage er- 
lebte. 

Die vorliegende Arbeit soll sich mit einem Teil- 
problem der geometrischen Raumakustik beschaf- 
tigen, das allerdings eine groBe Bedeutung hat, nim- 
lich mit der Verstarkung des direkten Schalles durch 
die ersten Rickwiirfe. Es ist seit langem bekannt, 
daB man die Lautstirke eines Schallvorganges am 
Ort des Zuhérers dadurch erheblich erhdhen kann, 
da8 man ihn in sehr kurzen Zeitabsténden ein oder 
mehrere Male wiederholt, was praktisch durch die 
Schallreflexionen an den Raumbegrenzungsflachen 
in der Nahe der Schaliquelle geschieht. Allerdings 
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miissen die reflektierten Schallanteile innerhalb von 
rund 50 ms nach dem direkten Schall ,,Schallquelle 
—Beobachter“ eintreffen, d.h. also die Wegdiffe- 
renzen zwischen direktem und reflektiertem Schall 
dirfen nicht mehr als héchstens 17m betragen. 
Diese Fragen sind kirzlich auf Veranlassung des 
einen Verfassers durch H. Haas [1] quantitativ 
studiert worden. Seine Ergebnisse beziehen sich auf 
die Versténdlichkeit von Sprache unter dem Ein- 
fluB eines Hinzelechos und bestatigen im wesent- 
lichen die oben genannten Zahlen, erweitern sie al- 
lerdings in dem von vornherein zu erwartenden 
Sinne, daf schwachere Reflexionen groBere Lauf- 
zeitdifferenzen zulassen, starkere Reflexionen da- 
gegen kleinere Laufzeitdifferenzen bedingen. Starke 
Absorption der reflektierenden Wand bei hohen 
Frequenzen erweitert gleichfalls den zulassigen 
Kchobereich zu groBeren Laufzeitwerten hin. DaB 
auBerdem die Sprechgeschwindigkeit in die Resul- 
tate eingeht, ist selbstverstandlich. Die Versuche 
von Haas wurden mit einer elektroakustischen An- 
lage gewonnen, bei welcher der direkte Schall und 
das Echo von je einem Lautsprecher kamen. Hine 
Magnetophonapparatur mit endloser Bandschleife 
diente zur Verzégerung des Echos gegentiber dem 
direkten Schall. AuBerdem konnten die Schallstarke, 
die Klangfarbe sowie die raumliche Richtung des 
Kchos in beliebigem Verhiltnis zu den Werten des 
direkten Schalles eingestellt werden. Inzwischen 
wurde diese Versuchsanordnung! so erweitert, daB 
an Stelle eines Echos mehrere Echos, und zwar bis 
zu neun, erzeugt werden konnen. Aus den noch im 
Gang befindlichen Versuchen mége bereits ein Er- 
gebnis erwéhnt werden, das fiir die vorliegenden 
Betrachtungen von Bedeutung ist. Kin Echo, das 
gleich stark wie der direkte Schall ist und eine Lauf- 
zeitdifferenz zwischen 50... 100 ms hat, stért sehr. 
Nach dem tibereinstimmenden Urteil verschiedener 
Beobachtergruppen kann diese Storung wesentlich 
gemildert werden, wenn zwischen direktem Schall 
und Kcho ein oder mehrere Echos eingeschoben wer- 
den. Man sieht daraus, daB eine Serie von Echos sich 
anders verhalt als ein Einzelecho und, was von 
vornherein wahrscheinlich ist, weniger stort. Auf 
der anderen Seite scheint aber auch das folgende 
Ergebnis festzustehen: Ersetzt man ein Einzelecho 
durch eine Anzahl sehr kurz aufeinander folgender, 
lautschwacher Echos, deren Gesamtenergie der des 
Einzelechos entspricht, so ist auch der subjektive 
Eindruck der gleiche. 

Die Riickwiirfe innerhalb der zulissigen Zeit er- 
hdhen aber nicht nur die empfundene Lautstarke, 


1 Bearbeiter: G. SCHODDER. 


sondern sie lésen auch eine ,,raumhafte’ Wirkung © 
beim Beobachter aus. Eine punktformige Schall- 
quelle scheint sich riumlich auszubreiten, ohne daB 
aber dabei, der Richtungseindruck verlorengeht. 
Gerade die eben genannten elektroakustischen Ver- — 
suche von Haas zeigten mit groBer Hindringlich- | 
keit, daB der Richtungseindruck durch den zuerst ° 
ankommenden Schall bewirkt wird und auch dann 
nicht verlorengeht, wenn das aus einer anderen | 
Richtung kommende, aber nicht spater als héch- 
stens 50 ms eintreffende Echo vielmals starker ist — 
(etwa bis zu 10 db). Dies gilt auch dann, wenn wie- 
derum das Einzelecho durch eine Serie kurzzeitig 
aufeinander folgender lautschwacher Echos ersetzt 
wird. Dabei scheint tibrigens ihre Gesamtenergie 
den eben genannten Wert von 10 db Pegelunter- 
schied noch etwas tberschreiten zu diirfen. 

Will man also mit Hilfe von zeitlich richtig auf- 
einander folgenden Riickwiirfen die Schallstarke des 
direkten Schalls erhohen, so mu8 man in bekannter 
Weise Wand und Decke in der Nahe der Schallquelle 
geeignet orientieren. Dies ist wohl zum ersten Male 
mit groBer Folgerichtigkeit im Salle Pleyel in Paris 
gemacht worden. Aber auch Hinzelreflektoren sind 
bekannt, z. B. die ebenen Kanzeldeckel in Kirchen. 
Auch je ein Reflektor[2] links und rechts von einem 
Sprecher oder ein groBer, parabel- und hyperbel- 
artig gekriimmter Hohlspiegel [3] hinter einem 
Sprecher sind gelegentlich versuchsweise verwendet 
worden. Hin Nachteil dieser Anordnungen ist, daB 
sie dem Auge des Beobachters unmittelbar ausge- 
setzt sind, ein gréBerer jedoch, daB man damit die 
seitlichen Wande des Raumes in der Nihe der 
Schallquelle nicht richtig ausnititzt. Es bereitet zwar 
keine Schwierigkeiten, den Seitenwainden im Grund- 
ri8 eine solche Stellung zu geben, da8 méglichst 
viel Schall nach den weiter riickwarts im Raum lie- 
genden Zuschauerplatzen reflektiert wird. Hine sol- 
che Anordnung wird bei vielen groBen Raumen in 
mannigfachen Formen durchgefihrt (Facherform, 
Birnenform, Trapezform usw. des Grundrisses). 
Diese Uberlegungen gelten aber nur fiir die Ebene, 
in der sich Schallquelle und Beobachter befinden, 
d. h. also im allgemeinen fiir eine horizontale Ebene. 
Der Schall jedoch, der auf die oberen Teile der Sei- 
tenflachen auffallt — und dies sind praktisch groBe 
Flachen, bzw. von der Schallquelle aus gesehen 
groBe Raumwinkel —, wird dagegen bisher nicht fiir 
eine erste Reflexion ausgeniitzt. Es liegt nun nahe, 
dies dadurch zu erreichen, daf man die betreffen- 
den Teile der Seitenflichen in eine groBe Reihe 
von ebenen Schallspiegeln auflést und diese um 
eine vertikale und eine horizontale Achse so dreht, 
da8 sie den Schall an den gewiinschten Ort, nim- 
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_ lich in die hinteren Raumteile dirigieren. Diese re- 


flektierten Schallanteile haben ebenso wie die an 
der Decke reflektierten Schallkomponenten den 
groBen Vorteil, daB sie nicht langere Strecken tiber 


das Publikum hinweglaufen und dadurch eine zu- 


sitzliche Schallschwachung bei mittleren Frequen- 
zen (etwa 0,5 db/m) erfahren. Ubrigens gibt es ein 
einfaches Verfahren, jeden dieser Schallspiegel zu 
justieren, sofern sie drehbar montiert sind. Man be- 
festigt auf den Platten einen groBen optischen Spie- 
gel und 1é8t das Licht eines Bihnenscheinwerfers, 
der in der Mitte der vorderen Bihne aufgestellt ist, 
daran reflektieren. Die Schallspiegel sind z. B. 
etwa 2 m? groB und kénnen aus gut reflektierenden 
Sperrholz- oder Bauplatten bestehen. Es kann aber 
auch angebracht sein, die Spiegel gleichzeitig als 
Schluckanordnung fiir tiefe Frequenzen auszubil- 
den. 

Vor die ganze Anordnung kommt ein optisch un- 
durchsichtiger, akustisch aber véllig durchlassiger 
Schirm, z. B. aus Drahtgaze, aus Rohrgeflecht, aus 
weitmaschigem Stoff, aus sehr stark perforierten 
diinnen Platten oder aus einem diimnen Holz- oder 
Metallgitter?. Damit ist jede Formgebung des Rau- 
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Um die Wegdifferenz zwischen direktem und re- 
flektiertem Schallstrahl klein zu halten, kommt es 
sehr darauf an, die der Schallquelle benachbarten 
Wandflachen so auszurichten, da8 ihre Spiegel- 
punkte in bezug auf die Schallquelle méglichst dicht 
beieinander liegen, und da} méglichst nur erste Re- 
flexionen auftreten. Abb. 1 stellt vier Beispiele zu- 
sammen, die der Anschaulichkeit wegen sich nur 
auf die Darstellung in der Ebene beschrinken; es 
sind dies der Fall der parallelen Seitenwainde mit 
einer ebenen Riickwand, der Fall der schragen Sei- 
tenwande mit einer ebenen Riickwand, der Fall ei- 
ner Parabel, sowie die Auflésung aller drei Wand- 
flachen in der Nahe der Schallquelle in eine Reihe 
von Hinzelreflektoren. Nur in den beiden ersten 
Fallen spielen Reflexionen zweiter oder hoherer 
Ordnung eine nennenswerte Rolle. Beim Hinsetzen 
der absoluten Zahlenwerte fiir die Raumdimen- 
sionen und bei der Annahme eines bestimmten Ortes 
fir den Beobachter darf der Abstand des Beobach- 
ters von jeder Spiegelschallquelle héchstens 17 m 
groBer sein als der Abstand des Beobachters von 
der Schallquelle*. Das gleiche gilt fiir den dreidimen- 
sionalen Fall, wobei, wie eben auseinandergesetzt 

wurde, die Auflosung gerade der oberen 

Teile der Seitenflachen in eine Reihe von 
Schallspiegeln groBe Bedeutung hat. 

Die statistische Theorie zeigt, da8 die 
Energiedichte Er im diffusen statio- 
naren Schallfeld von der Schalleistung 

a L der Schallquelle (Kugelschallquelle) 
und von der gesamten Schallschluck- 
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Abb. 1. Schallquelle und Bildpunkte der Schallquelle bei verschiedenen 
Raumbegrenzungen (parallele Seitenwande, schrage Seitenwande, 


Parabel und Spiegelanordnung). 


mes im Sinne des Architekten méglich. Ubrigens 
kann man an Stelle einer gleichmaBig geneigten und 
unmittelbar sichtbaren Decke auch hier dasselbe 
Prinzip der Unterteilung der Decke in viele, nach 
akustischen Gesichtspunkten geneigte Schallspiegel 
hinter einer schalldurchlissigen ebenen Abdeckung 
anwenden. 

2 Alle diese Anordnungen sind bis zu einer oberen Grenz- 
frequenz praktisch véllig schalldurchlassig. Man muB sie so 


bemessen, daB die Grenzfrequenz am oberen Ende des 
Ho6rbereichs liegt. 


flache A des Raumes abhangt, und zwar 
in der Form AL 
x aa Fae 
Andererseits ist die Energiedichte Ep 
des direkten Schalls 
15 Al 
Los 4n re" 
In graphischer Darstellung gibt dies 
Abb. 2a an, in der der Schallpegel uber 
dem Logarithmus der Entfernung r auf- 
getragen ist. 

Die im gr6éBten Teil des Raumes praktisch kon- 
stante Energiedichte Er des diffusen Schallfeldes 
kommt bekanntlich durch die zahlreichen Schall- 
reflexionen an allen Flichen des Raumes zustande. 
Mit Hilfe von Spiegelanordnungen gelingt es, den 
Horern zahlreiche niitzliche Schallreflexionen in- 
nerhalb einer Zeitdifferenz von 50 ms nach dem di- 
rekten Schall zuzufiihren. Die Energie, die in diesen 


3 Genaue Werte erhalt man entsprechend den Messungen 
von HAAS und den obigen Angaben iiber Mehrfachechos. 
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Reflexionen enthalten ist, ist im akustisch-physio- 
logischen Sinne ebenso wichtig wie der direkte 
Schall. Sie wird dem diffusen Schall entzogen, inso- 
fern sie vom Publikum absorbiert wird. Es seien « 


log E Ep 


emails ie log r 


Abb. 2. Schallenergiedichte des direkten und indirekten 
Schalles in Abhingigkeit von der Entfernung 
Schallquelle—Beobachter (ohne und mit Schall- 


spiegel); 
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der Bruchteil der von der Kugelschallquelle nach 
allen Seiten abgestrahlten Schalleistung L, der di- 
rekt oder innerhalb von 50 ms zum Horer gelangt, 
6 der Bruchteil des Raumwinkels 4x, unter dem das 
Publikum von der Schallquelle aus gesehen er- 
scheint, und ap die Schluckzahl des Publikums. 
Dann sind die durch die Spiegelanordnung ernied- 
rigte diffuse Energiedichte Ep, 
Eg= (aay) 

und die durch die Spiegelanordnung erhéhte Ener- 
giedichte Ep des direkten Schalles einschlieBlich 
seiner 50 ms-Trabanten 


, aL 1 
Samo 


Dabei ist r», ein mittlerer Wert der Abstiande zwi- 
schen Schallquelle bzw. Spiegelpunkt und Beob- 
achter. 
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Abb. 3. Schallriickwiirfe im 50 ms-Intervall und Nachhall 
bei Erregung des Raumes mit Impulsen. 


In Abb. 3 ist schematisch der zeitliche Schallver- 
lauf bei Erregung des Raumes mit Impulsen ange- 
geben, wobei die Riickwiirfe in dem Zeitintervall von 
50 ms nach dem direkten Schall einzeln eingezeich- 
net sind. Die nach dem kritischen Zeitintervall + 
erfolgenden Schallreflexionen sind in dem allgemei- 


nen Nachhall zusammengefaBt; ihre Hillkurve ist 
gestrichelt eingezeichnet. 

Je groBer man in den obigen Formeln «, d. h. also 
den Teil der Strahlung einer Kugelschallquelle 
macht, der das Publikum trifft, je mehr man also 
die Schallquelle richtet, um so mehr Energie ent- 
zieht man dem Nachhallschallfeld, und damit ni- 
hert man sich auch in geschlossenen Riéiumen dem 
Fall einer im Freien aufgestellten Orchestermuschel, 
der musikalisch keineswegs ideal ist. Man sieht an 
diesem Grenzfall, daB es auch méglich ist, die Ziele 
der geometrischen Akustik zu tibertreiben. Anderer- 
seits darf man vielleicht auch den SchluB ziehen, 
da8 in Réumen mit starker Betonung der kurz nach 
dem direkten Schall eintreffenden Reflexionen gré- 
Bere Nachhallzeiten als die bekannten, sogenannten 
optimalen Werte zulassig, eventuell sogar nétig sind. 

Zum SchluB seien einige praktische Anwendun- 
gen besprochen, die kirzlich ausgefiihrt worden 
sind. Beim Wiederaufbau der Hamburger Staats- 
oper* (Zuschauerraum 4400 m°, Biihnenraum 1600 
m%, 1100 Personen) mute aus auBeren Griinden 
ein Provisorium geschaffen werden, das aus der 
Anfiigung eines schmaleren zusatzlichen Zu- 
schauerraumes an den noch stehengebliebenen 
Bithnenraum bestand, der seinerseits zur Hilfte als 
Zuschauerraum, zur Halfte als Biihnenraum dient. 
Hierdurch ist die eigenartige Form des Grundrisses 
zu erkliren, wie sie aus den Abb. 4a und 4b hervor- 
geht. Die Dimensionen sind: Liinge 31 m, Breite im 
vorderen Teil des Zuschauerraumes 25 m, im riick- 
wirtigen Teil 16 m. 
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Abb. 4a. Grundri8 und Langsschnitt der Staatsoper Ham- 
burg. 


Die sonst immer vorhandenen Schallreflexionen 
an den Seitenwanden in der Ebene Bithne—Zuhérer 
fallen bei der vorliegenden GrundriBform praktisch 
aus, zumal an diesen Stellen die offenen Zuginge fiir 
das Publikum sowie die Garderobe sich befinden. 
Hier lag also ein Zwang vor, die oberen Flachen der 
Seitenwande fiir niitzliche Schallreflexionen zu ,,ak- 
tivieren“*. Dies geschah durch je 16 iiber die beiden 
vorderen Seitenflachen verteilte Spiegel von insge- 
samt 64 m? Fliche. Auch am vorderen Teil der Decke 


4 Architekt W. KALLMORGEN, 
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wurde ein groBer Spiegel 
angebracht. Diese genann- 
ten MaSnahmen haben 
sich sehr bewahrt. Zusam- 
men mit der Anwendung 
einer starken Uberhohung 
der Sitzreihen im Rang 
und im Parkett wurde 
eine ganz vorziigliche Ver- 
standlichkeit der Sprache 
auf allen Plaitzen erzielt, 
auch unter dem sehr weit 
vorspringenden Balkon. 
Im leeren Raum wurden 
_ Silbenverstandlichkeits- 

-Messungen durchgefihrt, 
die einen Abfall der Ver- 
standlichkeit auf den letz- 
ten Platzen von nur weni- 
gen Prozenten gegeniiber 
dem Wert in den ersten 
Parkettreihen ergaben. 


Abb. 4b, Innenansicht der Staatsoper Hamburg (Aufn. Wei’, Hamburg). 


mentes von Holz- und Metallplatten zeigten, daB 
| “te zwei diinne Aluminiumplatten von je 0,5 mm Starke, 

mit einer Pappezwischenlage zusammengeklebt, 
: einen ungefahr gleich groBen mechanischen Ver- 
~ | ohne Publikum lustwinkel wie Sperrholzplatten und den gleichen 
Spier ay an Schluckgrad haben. So wurden rund 400 m? dieser 
schwingungsfihigen Plattenabsorber in 5... 15 cm 
foe Abstand von der Raumwand mit einer Randdamp- 
i za r r fung des Luftpolsters (Auskleidung der Rander mit 


+ 
= 1100 Personen 
i} 


+ —— 


ema otha} einem porésen Schluckstoff) eingebaut. Die Nach- 
100 200 500 1000 2000 5000 10000 hallzeiten des Theaters ohne Publikum und bei 
Sa Titc) 


Nachhalizeit der Staatsoper 
Hamburg (mit und ohne Pu- 
blikum). 


Abb. 5. 


Der Vollstandigkeit hal- 
ber seien noch einige Mes- 
sungen in der Staatsoper 
erwihnt, die nicht in die 
geometrische, sondern in 
die statistische Raum- 
akustik gehéren. Aus 
feuerpolizeilichen Griin- 
den durfte in der Stuats- 

oper Hamburg die sonst 
viel angewendete Sperr- 
_ holzverkleidung nicht ge- 
| nommen werden. Einge- 

hende Versuche hinsicht- 


lich des Dampfungsdekre- Abb. 6. Mensasaal in Bonn (Aufn. Schafgans, Bonn). 
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Abb. 7a. Grundri8 des Landestheaters Hannover. 


voller Besetzung sind in Abb. 5 angegeben. Die 
akustische Anregung des Raumes erfolgte im leeren 
Theater mit Knallen, im besetzten Theater wurden 
die Messungen wahrend einer Opernauffiihrung bei 
Generalpausen nach Fortissimostellen des Orche- 
sters gemacht. Die Auf- 
zeichnung geschah iiber 
Oktavfilter mit dem Neu- 
mannschen Pegelschrei- 
ber®. Die Kurven in Abb. 5 
zeigen wbrigens, daB die 
Nachhallzeiten etwas iiber 
den Werten liegen, die 
man fiir die RaumgréBe 
der Hamburger Staats- 
oper sonst als optimal an- 
sieht, und daB sie nach tie- 
fen Frequenzen zu ziem- 
lich ansteigen. Die aku- 
stischen Eigenschaften des 
Raumes werden als gut 
bezeichnet. 


> Bei diesen Messungen und 
bei den Messungen der Sil- 
benverstandlichkeit wurden 
wir in dankenswerter Weise 
von der Akustischen Abtei- 
lung des Nordwestdeutschen 
Rundfunks in Hamburg un- 
terstutzt. 
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Ein zweites Beispiel fiir die Anwendung der Spie- 
gel ist der groBe Saal des Studentenhauses in Bonn 
(Abb. 6)&. Er dient sowohl als Restaurationsraum 
wie als Versammlungsraum und als Raum fir musi- 
kalische und kiinstlerische Darbietungen aller Art. 
Er ist rund 3000 m® groB und faBt 1000 Personen. 
Die raumakustische Schwierigkeit besteht darin, 
da er sehr lang (37 m) und schmal ist und wegen 
seiner Benutzung als Restaurant keine Uberhéhung 
der Sitzplatze zula8t. Um auch hier eine gute Ver- 
stindlichkeit zu erhalten, wurde vorgeschlagen, die 
frei verfiigbaren Flichen in der Nihe der Schall- 
quelle fiir eine Spiegelanordnung auszuniitzen, und 
zwar sind dies die Riickwand und die Decke. Bei 
der Riickwand ist die Spiegelanordnung nicht ver- 
kleidet, sondern frei sichtbar, wahrend die Decke 
auf ihrer ganzen Breite und in etwa 10m Lange, 
von der Riickwand aus gerechnet, mit Spiegeln ver- 
sehen ist, die um die horizontale und um die verti- 
kale Achse so gedreht sind, daB modglichst viel 
Schall auf die riickwartigen Plitze im Parkett und 
auf die Galerie kommt. Die Deckenspiegel sind 
durch ein Drahtgewebe verkleidet. Der Erfolg dieser 
Ma8nahmen ist gut. Auch fiir die letzten Reihen des 
Saals wird die Sprechverstindlichkeit als gut be- 
zeichnet. Eine vorhandene Lautsprecheranlage 
wird nur bei sehr leisen Sprechern eingeschaltet. 


6 Architekt T. KLEEFISCH, Bonn. 


Abb. 7b. Innenansicht des Landestheaters Hannover 


(Aufnahme Bauverwaltung, Hannover). 
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_ Zum SchluB sei noch das im Ausbau befindliche 
Landestheater Hannover’ erwaéhnt (Abb. 7a und 
Tb). Der Zuschauerraum hat ein Volumen von 
7000 m3. Infolge des hier wohl zum ersten Male in 
die Theatertechnik eingefiihrten Rundvorhanges 
und infolge der auf der vorderen Biihne aufgestellten 
-zahlreichen Beleuchtungsscheinwerfer fallen gegen- 
iiber der sonst iiblichen Theaterbauweise (Prosze- 
niumswande) erhebliche Teile der unteren Seiten- 
_wiinde fir eine Schallreflexion aus. Es sind dies aber 
gerade die Seitenwande in der Nahe der Schall- 
quelle, die besonders wichtig sind. Um diesen Man- 

gel bis zu einem gewissen Grade auszugleichen, 
wurden auf der Vorbithne an der in Abb. 7 bezeich- 
-neten Stelle zwei Reihen von Schallspiegeln, die 
sich vom BihnenfuBboden bis zur Decke erstrecken, 


7 Architekten W. KALLMORGEN und Cl. HOFFMANN. 


angebracht. Auch hier kann man bei Personen, die 
auf dem vorderen Teil der Buhne sprechen, sehr gut 
den Einflu8 der Spiegel in einer Zunahme der Laut- 
starke und der Verstindlichkeit feststellen. Leider 
muBBten aus anderen Griinden die Spiegel schrag 
aufgestellt werden, so daB sie fiir Schauspieler, die 
weiter hinten auf der Bihne stehen, nicht mehr 
wirksam sind. Eine zur Bihnenachse parallele Auf- 
stellung der Spiegel wiirde ihren Hinflu8 erheblich 


erhohen. 
(Eingegangen am 15. Juni 1951.) 
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DESIGN AND ANALYSIS OF SUBJECTIVE ACOUSTICAL 
EXPERIMENTS WHICH INVOLVE A QUANTAL RESPONSE 


by D. L. RICHARDS, B. Sc., F.S.S., A.M.I.E.E. 
Post Office Research Station, Dollis Hill, London 


Summary 


Subjective acoustical measurements involve quantities which do not possess an intrinsic scale 
of magnitude and so the only observations possible are in the form of very simple classification 
judgments or quantal responses. These measurements are usually also concerned with relating the 
subjective response to some objective, physically measurable stimulus. Much of the methodology 
developed in the broader field of biological experimentation is surprisingly appropriate to these 
problems. Some examples are given of the application of ‘‘probit analysis” and ‘‘analysis of vari- 
ance” to the measurement of telephone circuits. These statistical methods used with appropriate 
designs of experiments not only yield valid confidence intervals but also point to the most efficient 
method of conducting the measurements through separation of the several sources of variation. 


Sommaire 


Les mesures acoustiques subjectives font intervenir des grandeurs dont les amplitudes ne 
peuvent pas étre rapportées & une échelle intrinséque; leg observations faites ne peuvent étre 
exprimées que sous la forme d’une trés simple classification («quantal responses»), Les mesures 
ont généralement trait, par ailleurs, & la relation existant entre une réponse subjective et un 
certain stimulus objectif et physiquement mesurable. Une grande partie de la méthodologie mise 
au point pour important domaine de l’expérimentation biologique convient étonnamment bien 
a l’étude des problémes ci-dessus. On donne ainsi quelques applications de la «probit analysis» 
et de l’<analysis of variance» 4 1’étude des circuits téléphoniques. De telles méthodes statistiques, 
appliquées & des schémas expérimentaux appropriés, non seulement permettent d’obtenir des 
inégalités sires, mais aussi conduisent au procédé de mesure le plus efficace, en ce qu’elles per-. 
mettent de séparer les différentes causes de variation. 


Zusammenfassung 


Es wird gezeigt, daB gewisse statistische Methoden, die in der Biologie verwandt werden, sich gut 
fiir subjektive akustische Messungen eignen. Einige Beispiele fiir die Anwendbarkeit dieser 
Methode auf dem Gebiet des Fernsprechwesens werden angefiihrt. 


however, has an impetus been given by the devel- 
opment of appropriate methodology in the much 
broader field of biological experimentation. These 
methods are now being used in telephone meas- 
urements. As they are also applicable to a very 


1. Definition of the class of problems covered 


The technique of assessing properties of telephone 
circuits by subjective responses has been gradually 
developing over several decades, Only recently, 


| 
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wide range of subjective experiments, the problems 
are treated here in a rather general manner. 

Subjective acoustical measurements differ from 
most ordinary physical measurements in an impor- 
tant respect, namely that the quantities which are 
being measured do not possess an intrinsic scale of 
magnitude. For example, whereas there is no diffi- 
culty involved in the conception: “‘Jength A is 1.4 
times length B’’, the phrase: ‘‘telephone circuit A is 
1.4 times as easily usable as telephone circuit B’’ is 
virtually meaningless. If a subject is presented with 
two circuits, the only relevant subjective observa- 
tion that it is possible to make is to ask him to de- 
cide upon one of a number of possible responses to a 
question such as, for example: “Which of the fol- 
lowing responses is the most appropriate ?”’ 

(1) A is more easily usable than B, 

(2) Bis more easily usable than A, 

(3) A and B are equally easily usable. 

It is usually possible to restrict the choice to one 
of two responses, either by combining (3) with (1) or 
(2) or by disallowing (3). The methods described be- 
low are appropriate when there are two alternative 
responses. Such “‘yes/no”, “‘all/none’’, “‘positive/ 
negative” or “death/survival’’ responses are termed 
“quantal” responses. 

If three possible responses must be allowed the 
method to be described can still be used but the 
analysis becomes more complicated. The results 
would then be analysed (a) with (3) combined with 
(1) and (b) with (3) combined with (2); the interval 
between the results for (a) and (b) then gives a 
“range of indifference’. 

As is shown later, quantal responses are also in- 
volved when measurements are conducted by mak- 
ing some adjustment so that, for example, equal 
loudnesses are obtained. 

A second feature of subjective acoustical meas- 
urements is that one is usually concerned with re- 
lating the subjective response to some objective, 
physically measurable stimulus. The level of such 
a stimulus is then used as a measure of the corre- 
sponding subjective response. 

Two kinds of experiment will be distinguished 
here: the first uses the judgments of a large number 
of subjects, each subject being used once only; the 
second makes use of a crew containing a limited 
number of subjects, each member of the crew being 
used many times. There exist different situations 
where either one of these is more appropriate. For 
example, if completely unpracticed subjects must 
be employed the first kind must be used; if, however 
great precision is required it would be better to use 
the second kind. 


2. An experiment using a large number of subjects | 
each once only 


A telephone assessment problem will be taken as 
an example of the use of a large number of subjects. 
In this case generality of the results was more im- | 
portant than precision. The problem was to deter- 
mine the minimum received level of a telephone 
circuit which would be regarded as usable by at 
least 50% of subjects. The subjects were asked to 
converse in pairs over the circuit set to a given level 
and then to give their decision by answering the 
question: “Was the circuit usable or not usable?” 
Seven different values of level were used altogether 
in the whole experiment, but each pair of subjects — 
was given one value only. The subjects of each pair 
gave their decisions independently by signal with — 
a push-button. A total of about 850 subjects was 
used and the percentage of positive (usable) respon- 
ses determined for each level. The level, x, was ex- | 
pressed in logarithmic units (decibels). The results 
which apply for the subjects who used one particu-— 
lar end of the circuit and who were referred to as 
“callers” are given in Fig. 1. The results for the sub- | 
jects who used the other end were treated separately 
lest the two sets were not independent. 
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Fig. 1. Percentage positive response (usability) as a func- 
tion of sensitivity. 


2. 1. Analysis of the results 


The results given in Fig. 1 will now be considered 
in detail as an example of the application of statisti- 
cal methods which have been developed in the bio- 
logical field for assaying the potency of drugs. Such 
work is termed “biological assay” and the partic- 
ular statistical method used here is termed “probit 
analysis’. 

When analysing any experimentally obtained 
data it is necessary to make certain assumptions [1]. 
These are acceptable either because they are regard- 
ed as self-evident or in view of previous experience. 
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The validity of any method of analysis depends on 


the truth of its assumptions and, in general, if any 


of these are false, conclusions drawn from the data 
will be unfounded and may be misleading. It is for- 
tunate, however, that some of the necessary assump- 
tions need only be “approximately”’ true; that is to 
say, if the data themselves do not render them un- 
likely they may be accepted. If this were not the 
case scarcely any statistical methods could ever be 
used in practice! 

The first assumption which concerns the methods 
to be described here is that, if the number of ob- 
servations were increased indefinitely and the inter- 
vals between the values of x made smaller (while 
still covering the same range), the observed points 
would lie on a smooth curve. This may be regarded 
as self-evident. The second is that the observations 
actually made form a random sample of the infinite 
population of possible observations which might 
have been made. This may, under certain conditions 
defined later, also be accepted. The third assump- 
tion is that the smooth curve referred to above has 
a particular mathematical form. 

In biological assay terminology, a stimulus is 
applied to the subject and the response is measured. 
Tn the case of quantal responses, if P is the propor- 
tion of positive responses and «x is the stimulus, P is 
a function of x. Many mathematical functions have 
been suggested [2], [3], but only the probability 
integral will be considered here. (The word probit is 


_ derived from this function.) There is no special rea- 


son for using it here except that the necessary tables 
are available and that the data are not in disagree- 
ment with it. It also so happens that practically 
identical results are obtained if one of the alterna- 
tive functions is used. This particularly applies if 
only the value of x is required which corresponds to 
P = 0.5. In probit analysis it is assumed that: 


P= : 2 eh dl ERE 
y2rr 


where p. and 2 are to be estimated from the data. 
The assumption of a function P= f(x) can be con- 
sidered from a different aspect as follows. Let us 
consider a hypothetical experiment using a given 
subject; the value of x is gradually increased until 
the subject just gives a positive response. This criti- 
cal value is termed the “tolerance”. If we now 
assume that the tolerances for the population of 
subjects are distributed in the form df(x)/dx, this 
is identical to the former assumption. If f(x) is the 
probability integral, then, of course, the distribution 
is Gaussian. In almost all statistical analysis of ex- 
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perimental data a Gaussian distribution is assumed 
without any special positive evidence in its favour. 
As will be shown later the conception of tolerances 
is very convenient when complicated experiments 
are considered. 

The analysis therefore consists of choosing the 
“most likely’? values for the two parameters of the 
function. This estimation is carried out by the meth- 
od of maximum likelihood. Details of the compu- 
tational procedure are given elsewhere [4] but a 
sketch of the principles will be given here. It is con- 
venient to write the function relating the proportion 
of positive responses as P=f(x; A,B) where A and 
B are two parameters which define p and 4. 

When x=~,, let P=P, and let n, be the number. 
of positive responses observed out of N,. Then the 
probability that exactly n, positive responses would 
be obtained at random out of N, trials is 


n,!(N, —n,)! 


The compound probability that exactly n, would 
be obtained out of N,, that n, would be obtained 
out of N,, and so on, is the product of all the similar 
terms. This product (or strictly a quantity propor- 
tional to it) is termed the “likelihood” of the set of 
observations. The method of maximum likelihood 
consists of so choosing the parameters that this 
quantity is a maximum. The calculations yield m, 
an estimate of u, and b, an estimate of 1/4. mis the 
median value of x and 2 is the standard deviation 
of the tolerance distribution. It is also possible to 
determine confidence limits for m and 6. A com- 
plete fitted curve may be drawn to the observed 
points giving the value of x for any desired value of 
P. Such a curve has been drawn on Fig. 1 and its 
parameters are m= 25.8, 4=12.1 (both in decibels). 


P, eu py a) ; 


2. 2. Assumptions implied by the method of analysis 


In addition to the assumption of a particular 
mathematical form for P as a function of x, other 
assumptions are made implicitly when the probit 
method of analysis is used. The most important of 
these is that the reactions of the subjects are inde- 
pendent of the batches into which they are divided 
for allocation to the various levels of x. To eliminate 
possible sources of bias it is necessary to take certain 
precautions such as: to ensure that the subjects do 
not know the value of x which is being used; to allo- 
cate the subjects at random to the values of x. The 
importance of randomization cannot be too strongly 
emphasized. If correct randomization is not carried 
out the results may be seriously biased and wrong 
conclusions may be drawn. 
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The dispersion of the observed results about the 
fitted curve can be tested statistically and this test 
should normally be included in the computations. A 
statistic y? is calculated and, if 7? is not significant, 
the dispersion may reasonably be attributed to 
chance. If, however, 2 is found to be significant, 
either the probability integral was not a suitable 
choice of function or (more frequently) the reactions 
of the subjects were not independent of x. The 
latter effect may arise through a systematic, though 
inadvertent, allocation of subjects. For example, 
the values of x may have been used in a systematic 
order which, perhaps unknown to the experimenter, 
happened to coincide with the order of some other 
factor. It may have been that there were, say, seven 
values of x and perhaps the subjects were called in 
groups of seven so that the subject who had to wait 
the longest always used the same value of x. This 
is a fairly obvious fault of organization which. would 
usually be avoided in any case, but it is impossible 
to foresee every factor which may conceivably af- 
fect the results. The difficulty can be completely 
overcome, however, by using the values of x in 
random order. 

It is sometimes objected that the use of randomi- 
zation increases the dispersion of the experimental 
results. This is sometimes true but the reduced dis- 
persion by not randomizing is then spurious and 
will give a false suggestion of precision unless the 
systematic arrangement is in accordance with a pro- 
perly planned, executed and analysed experiment. 
An example of this type of experiment is given 
later. 


3. Measurements using the same subject repeatedly 


It is frequently necessary in subjective experi- 
mentation to use a comparatively small team of 
subjects each member of which is used for every 
condition and, in order to increase precision, the 
entire set of measurements is repeated several times. 
This procedure requires a very careful consideration ; 
there are circumstances under which increasing the 
number of observations does not increase the pre- 
cision as much as might be expected. In such cases 
inappropriate analysis can lead to misleading con- 
clusions. 

As an example a very common problem in tele- 
phone assessment will be considered. There are 
two telephone channels and one of them is adjust- 
able in sensitivity. It is required to determine the 
sensitivity setting which causes the loudness of the 
received speech over both channels to be the same. 
Usually one of the channels is a Reference Channel 
(R) and the other is to be assessed in terms of the 


Reference Channel. Sometimes the channel to be 
assessed (X) is the adjustable one and the setting of 
X is required which corresponds to the fixed Refe- 
rence Channel. In other instances the Reference 
Channel is adjustable and its setting is required — 
which corresponds to the fixed channel X. The 
latter situation is analogous to audiometry and — 
subjective loudness measurements. 


3. 1. Three methods of loudness balancing 


The method of loudness balancing to be described 
is novel not only as regards the way in which the 
data are analysed but also the actual conduct of the 
observations. 

There are two commonly used methods which 
will be briefly described first. Method I (see Fig. 2) 
uses a team of three persons; on acts as talker, one 
as listener and one as balancer. The talker provides - 
a constant level of speech by repeating a phrase 
such as “ONE-TWO-THREE-FOUR-FIVE”. The — 
listener is required to signal his decisions to the 
balancer; each of these decisions is either “X is 
louder than R” or “R is louder than X”. The | 
decision relates to the loudness of the speech 
received through the two channels. The task of 
the balancer is to adjust the sensitivity of, say, X 
(if X is the adjustable channel and R is fixed) in 
accordance with the decisions of the listener until 
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Fig. 2. Three methods of loudness balancing. 
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he has decided what setting of X corresponds to 


equality of loudness as judged by the listener. The 


task of the balancer is very difficult because he is 
required to have insight into the mind of the liste- 
ner. Considerable training is necessary when this 
method is used. 

Some reduction of the amount of skill required is 
achieved by using method II (see Fig. 2) where the 
tasks of listener and balancer are combined. Only 
two persons are necessary, a talker and a listener. 
The talker acts as in method I but the listener is 
given control of the sensitivity of X which he must 
adjust until he has decided that the speech from 
both channels is equally loud. The essential part of 


any subjective measurement is the making of de- 


cisions. Both the methods so far described require 
also some action consequent upon each decision. 
This latter element is undesirable and can with 
advantage be eliminated so that all action such as 
the setting of the sensitivity control is determined 
in advance. Method III limits the subjective ele- 
ment of the measurement to the mere giving of 
decisions. 

Method III (see Fig. 2) requires two persons 


| whose tasks are as follows. The talker, as before, 
_ provides the speech input to the telephone channel, 
_ but he also sets the sensitivity to a predetermined 


value. The listener gives his decision which is record- 
ed. The talker then sets the sensitivity to the next 
value and the listener again gives his decision. The 
observations are therefore obtained as a set of sen- 
sitivity settings with their corresponding decisions. 
The listener is required to give about eight decisions 
in succession, the sensitivity settings being present- 
ed in random order. It is convenient still to use a 


team of three persons so as to have six different 
_ talker/listener combinations available. In a typical 


experiment eight different sensitivity settings might 
be used and each talker/listener combination be 
used twice; there would then be twelve decisions 
per sensitivity setting. Such a set of observations 
will usually suffice since it will yield a 95% confi- 
dence interval of about +- 0.8 db. The results ob- 
tained by method III are analysed by the probit 
technique and so attention must be paid to the con- 
ditions for validity. The most important require- 
ment is that the circuits and sensitivity settings 
must be allocated at random to the talker/listener 
pairs. A satisfactory method of doing this randomi- 
zation is as follows. Suppose that there are two 
telephone circuits, R and X, and that it is desired 
to adjust the sensitivity setting of X so that the 
loudness of the received speech is the same for each. 
There are eight settings of the sensitivity control 
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and six talker/listener combinations are to be used 
twice each, so that there will be twelve sets of de- 
cisions. The twelve sets are numbered and one of 
these numbers is drawn at random, with replace- 
ment, for each of the ninety-six decisions (twelve 
for each of the eight sensitivity settings). This re- 
sults in unequal numbers of decisions to be made by 
each talker/listener pair on each occasion. There is 
no disadvantage in this; it may even be an advan- 
tage that the listener does not know which decision 
is to be the last of a set. The results of such an ex- 
periment are given in the following Table I. In this 
experiment a given circuit (X) was adjustable and 
the setting of this adjustment was required which 
corresponded to a level of loudness of the received 
speech defined by the fixed setting of a reference 
circuit (R). The talker repeated a phrase of meaning- 
less syllables. The sensitivity settings of X are at 2 
db intervals, increasing in attenuation from | to 8. 
In the table the setting number is given together 
with the decision which was given; R indicates that 
the reference circuit was judged to give the louder 
received speech and X, the circuit being adjusted. 


Table I 
talker/ umber of sensitivity settin 
oior | Ustener arta Hegel 
decisions | nation decision given 
—— 
1 | A/B | 4X 4X 4X 6R 6R 5X 
2 | B/C | 4X 1X 4R 7R 4X 8R 
3).))) C/A) 2X 8X 3X TR SR 1X 3X.8X 
4 | A/C | 1X 7R 2X 7R 2X 8R 3X 6R 8R 
5 | C/B | 3X 2X 8R 6R 2X 7R 1X 2X 
6 | B/A | 6R 8R 1X 7R 3X 3X 2X 
7 | A/B | 7R 4R 6R 8R 8R 3X 6X 2X 7R 
8 B/C | 1X 2X 1X 6R 5R 7R 3X 2X 
1X 4X 56R 3R 
9 C/A | 6R 2X 6R 8R 7R 
10 A/C | 6R 8R 5R 4X 3X 5R 8R 
11 C/B | 5R 6R 5X 4X 6R 8R 6X 1X 4K 
12 B/A | 7R 1X 1X 4X 3X 6X 5R 1X 2X 7R 


These results are analysed by the probit method 
as before. An R decision is counted as a positive 
response and when the numbers of these are classi- 
fied under the eight sensitivity settings Table IT is 
obtained. 

Probit analysis of these data yields the following 
results: 
Median sensitivity setting (corresponding to 
50% positive response) =m=70.47 db. 
95% confidence interval of m=-+ 0.81 db. 
A = 1.94 db. 
The goodness of fit of the function whose par- 
ameters are m and 2 is satisfactory as judged by 
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y?. That is the departures of the observations 
from the fitted curve are within those to be 
expected by chance. 

The statement of a confidence limit always re- 
quires some explanation of its meaning; in the 
above it applies to repetitions of the measurements 
using the same talker/listener combinations allo- 
cated at random in the same way. 


Table II 
sensitivity total num- | number 
setting ber of of positive Porcen tee 
(attenuation) decisons | responses DEP vaEryS 
Ho. rrre | N a responses 
Lele y68 12 Oia bs 0 
2 GB a saliva ies ity nO) back ia 
3 67 12 1 8 
4 ain a 2 17 
5 | 71 12 7 58 
6 73 12 11 | 92 
7 15 12 12) se 100 
Ag Pei 12 12, 100 


3. 2. Effect of the subject 


The method given above satisfies the requirements 
for validity of the probit method and provides a 
test of the departure from a Gaussian distribution 
of the tolerances. The confidence interval can only, 
however, be calculated which applies to repetitions 
of the measurements with the same combinations of 
talker and listener. If the more general confidence 
interval is required which applies for any talker/ 
listener combinations drawn from the same popula- 
tion, a different experimental design must be used 
because the effect of talker/listener combinations 
must be estimated. It may, furthermore, be suspect- 
ed that there is an interaction of circuits and talker/ 
listener combinations, i.e. that 
a different result may be 
obtained for the difference 
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measurements must be so conducted that a separate 
estimate of m is obtainable from every set of (say, 
eight) decisions. These values of m form the units 
which are subjected to analysis of variance. | 
An estimate of m can be made from a set of, say, 
eight decisions given for eight different, equally 
spaced, settings of sensitivity. For example suppose 
the following results were obtained for one set of 
decisions: , 


setting [db] 0 2 4 6 8 101214 
decision XX XRRBRR 


A good estimate of m is the mean of the highest 
sensitivity which yields a negative response and the 
lowest sensitivity which yields a positive response 
[4], [5]. In this case the mean is of 4 and 6, namely 5. 
Sometimes there is no setting of sensitivity below 
which only negative responses have been given and 
above which only positive responses have been 
given. In such cases there may be said to be a cross- 
over. An example is as follows: 


setting [db] 0 2 4 6 8 10 12 14 
decision xX X RXR ROR 


The same rule may still be applied; a good esti- 
mate of m is the mean of 6 and 4, again 5. Actually 
the probability of occurrence of a cross-over is a 
function of the ratio d/’, where d is the interval 
between sensitivity settings. In this case d = 2 db. 
It is shown later that the proportion of cross-overs 
can be used to estimate 4. ; 

An example will now be given to illustrate this — 
method. Table III shows part of the results (those 
for two circuits only) of an experiment to determine ~ 
the attenuation setting of a reference circuit, R, 
which corresponds for loudness of received speech 
to each of five circuits, X,, X,, X;, X,, X;. The indi- 


Table III 


between two circuit pairs talker/ shee Se chores 
when two different talker /list- block syayns attenuation of R [db] attenuation of R [db] 
These questions can best be TA ol ty ee Peat op RRRERE RX oe 
answered by arranging the ex- | B/C RRAXXXXXX | BRE Xo 
periment so that the results C/A RRRRRRRR | RRRRXRXX 
can be analysed by the analysis a : = e m = - co = 3 = = = = se 
of variance method. It has B/A RRERERRXXX |RRRRRRRX 
been found that one talker/ pal aca 
listener combination which has m| Ap RRRRXXxx RRRRRRXX 
given one consistent result on B/C BBG Re se eee RR XxX Xz Sees 
one occasion will sometimes QA | RXXXXXXKXK | RRRXRXXX 
give a different, but still con- an Bo aR SS ne 
sistent, result on another oc- B/A RRRXXXXX |RRRXXXXX 


casion. For this reason the 


vidual values of m estimated as above have been 
entered in Table IV which also shows the summa- 
tions necessary for computing the sums of squares 
required in carrying out an analysis of variance. 
The observations were made in the following way: 
(a) The observations of Block I were completed 
before Block IL was commenced. 
(b) The talker/listener combinations were used in 
rotation in the order given. 
' (c) The circuits were used in random sequence. 
(d) The attenuator settings within each set of 
eight decisions were presented in random 
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There are two entries in Table IV, marked * and 
t; in the case of * all eight decisions were R (this 
appears on Table III); the mean attenuations have 
here been estimated on the assumption that the de- 
cision for a setting 2 db higher than the highest ac- 
tually used would have yielded a decision X. The 
entry marked + corresponds to a missing observa- 
tion; the value given has been fitted by “missing 
plot” technique. 

The results of the analysis of variance carried out 
on data given in Table IV are summarized in Table 
V. The computations for such a table are described 


order. elsewhere [6]. 
The following conclusions can be drawn 
Table IV from this analysis: 
; (a) there is no doubt that the circuits 
talker/ ancats ae pa yield truly different results, 
block | listener, 36 Apea ee) ctrnuita (b) the results with different talker/list- 
era 1 2 3 ire etaeee ener pairs are genuinely different, 
I AB 19 95 42 47 56 189 but these differences are not the 
B/C ll 19 39 4] 52 155 same in the two blocks; so much so 
C/A 23* | 25 42 45 60 195 that the differences between talker/ 
A/C WD] 19 | 8204 485) 54 158 listener pairs lack all consistency. 
ed coool Si pdo de wien) “Ags This is illustrated by the “very highly 
| ahs, ee ee Se significant” talker/listener pairs and 
ate oNEGUE ES) EE 1078 blocks interaction (denoted by***), 
II A/B 15 27 40 45 56 183 coupled with the non-significant 
B/C M5 hO 1. 84 a Be ase talker/listener pairs effect (tested 
\ | es | i 3 ro | 45 P | ie against the interaction); 
CB 13 27 40 | 49 | 56 185 (c) the components of variance enable 
I eins HAS walle al 3 | 21 34 | 41. | 54 163 confidence intervals to be calculated. 
Eau i7en | ado Mase "| 264 (308° 10a” The calculation of confidence intervals 
ean Sf requires some explanation since different 
Te = " oe ie | ih . : ue . ate values apply according to the purpose for 
hiscks: VL oya 39 48 89 sg |i | 362 which they are required. In principle four 
AUKeea 622 | + 42 70 88 | 110 332 values can apply in the present example 
c/B |, 30 OF fa, 0S 96 | 114 372 and these can be classified as follows: 
| B/A | 30 mO.. |, 88 92 | as 348 bis 1. for a circuit mean as an absolute 
sum | 174 | 284 | 446 | 538 | 662 | 2104 quantity, 
2. for a circuit mean relative 
to the mean of the group of 
Table V circuit means, 
| ke | | vasinnas:, | component. together with further classi- 
factor eh oe eee Meir | ais oe fication according to whether 
freedom | squares _ Sie ea F __ variance the confidence intervals are to 
sR AEE REE ERO EEE eEEEn EE ETEE EE apply for 
circuits (Q) 4 | 12681 | 3170.25 ST9*e* ere a aenoniaernenin of 
talker/listener : . 
pairs (C) | eo 246 49.20 | 8.99 | ViC)=0 talker/listener pairs used 
Blocks 3 tow AB) nbs 1 th gh, | one | EES” or(b) a group of talker/listener 
cael he | | pairs drawn at random 
‘ (QC) 20 105 | 5.25 | 0.96 |: VIQC) = 0 a 
i (QB) 4 35 | 8.75 | 1.60 ViQB) =0 om the same popula- 
(CB) 5 136 27.20 4.97*** | V(CB) = 4.3 tion of talker /listener 
Error 19 104 5.47 ) Vy = 5.5 pairs. 


¢ 


9) 


Table VI shows how these four values are calcu- 
lated. The table gives the coefficients of the terms 
which must be added together to give the variance 
of a condition mean. The standard error is the square 
root of the variance and the 95% confidence interval 
is -|- 1.96 times the standard error. In this table 7,. 


denotes the number of talker/listener pairs (in this 
case 6) and n,, the number of blocks (in this case 2). 


Table VI 
type of V vic) Wo Vv 
confidence. —° — M(GBy V(B) HK 2) (QB) i. (QC) 
interval "CNB Noha, NB ne NR no 
Ja ] 1 1 0 1 0 
tb l l 1 1 1 i 
2a, 1 0 0 0 1 0 
2 1 0 0 0 1 f 


The analysis of variance of the present set of re- 
sults gives V) = 5.5, V(CB)=43 and V(B) = 
V(C) == V(QB) = V(QC) = 0. The confidence inter- 
val of a circuit mean as an absolute quantity (la or 
Ib) is therefore -|- 1.96 y'5.5/12 +4.3/12 = +1.8 db. 
The corresponding value for a circuit mean relative 
to the group of circuit means is + 1.96 y 5.5/12 
— -|- 1.3 db. If two of the circuit means are to be 
compared with each other the confidence interval 
ot the difference between two circuit means is | 2 
times the latter confidence interval, namely + 1.9 db. 

The above example is of a comparatively small 
set of observations where only three factors are in- 
volved, (Q), (C) and (B), together with their inter- 
actions. lt is frequently necessary to include other 
factors such as microphones and receivers; the de- 
sign of the experiment then becomes correspond- 
ingly more elaborate. A suitable design, however. 
will enable all the appropriate confidence intervals 
to be estimated and, in addition, will provide a con- 
trol over the conduct of the measurements. For 
example, in the above a significant circuits and 
blocks interaction would have cast doubts on the 
stability of the physical circuits or on the conduct 
of the measurements. (For example a wrong circuit 
might have been set up on one or more occasions.) 

In Table III there are two sets of decisions 
which exhibit a cross-over (see section 3.2 above); 
out of the whole experiment containing sixty sets 
of decisions there are five cross-overs, i. e. the pro- 
portion of cross-overs is 0.083. The proportion of 
cross-overs is in a sufficiently large sample a func- 
tion of the ratio d, the sensitivity interval, to 4, the 
standard deviation of the tolerance distribution. 
(The actual proportion observed in any small sample 
is, like any other statistic, subject to sampling 
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variation.) In the present instance the value of d/r 
which corresponds to a proportion of 0.083 is 1.4; dis 
2dbso that 2 must be 2/1.4= 1.4 db. The variance of 
individual values of m estimated from sets of de- 
cisions as in this experiment is 0.57 Ad in this in- 
stance 1.6. This variance is significantly smaller 
than the error variance estimated from the analysis 
of variance, namely 5:5. The conclusion which may 
be drawn is that the variability of each subject is 
much less within a set of decisions than it is bet- 
ween sets of decisions. It follows that there is a limit 
to the precision with which it is useful to measure 
the individual values of m; indeed in this experi- 
ment four decisions at 4 db intervals (d = 4dh) 
would have resulted in an error variance of no more 
than 7.1 with a saving of half the number of deci- 
sions!; the precision would have been 77% with 
50% of the expenditure of decisions, an increase of 
efficiency of 54°%,. This conclusion is rather unex- 
pected but it illustrates the advantage of proper 
design. conduct and analysis of the measurements. 
The use of methods I and IT of conducting the loud- 
ness balancing would never have shown up this 
effect and much effort might have been employed in 
a futile attempt to improve the precision of the 
means for the circuits by conducting the individual 
balances more precisely. 

Table VII is given here which shows the rela- 
tionship between the expected proportion of cross- 
overs for values of d/i from 0.5 to 2.4. Certain ap- 
proximations have been made in the caleulations 
but the accuracy is quite adequate for the purpose 
of estimating 4. 


Table VII 


P is the proportion of cross-overs in a large sample 


d | Pp d P d Pi d p 
oe A A A 

05 0.540 1.0 0.196 15 , 0.070 le 20g Geer 
0.6 0,448 1.1 0.159 1.6 |.0.057! 2.1 7; 0.020 
0.7° 10.366; 1.2. 0.129, “1:7 C0046) 2s 
OS 10.298" MS FP O.205 1.8 | 0.038; 2.3 | 0.014 
0.9 ' 0.242. 1.4 0.086 1.9...:.0.03L. 224 


0.011 


4. Number of subjects or size of crew required 


When each subject is used once only or when 
measurements are so conducted as to yield results 
in the form of Table I, it is not difficult, given the 


1 d would then exceed yerror variance and so it would 
be advisable to arrange that there are still eight possible 
settings available at 2 db intervals of which only four 
alternate settings are used for any one set of decisions 
Nos. 1, 3, 5 and 7 (in ascending order of sensitivity) or 
Nos. 2, 4, 6 and 8 should be chosen at random. 
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value of 4, to calculate how many observations 
must be made to ensure that the confidence interval 
of each condition mean is reduced to a given value. 
When, however, the measurements are conducted 
in the form of a designed experiment with a syste- 
matic arrangement of such factors as talker/listener 
pairs and blocks, determination of the necessary 
number of observations depends not only on A but 
also on the magnitude of each component of vari- 
ance. For example, suppose that V(C), the talker/ 
listener pairs variance, only is significant and that 
the relevant 95% confidence interval is that which 
applies for a circuit mean as an absolute quantity 
and for a group of talker/listener pairs drawn at 
yandom from the population of talker/listener pairs 
(type 1b of Table VI). This confidence interval is 
-+ 1.96 V Vi/ncnz + V(C)/nc, where nc is the num- 
ber of talker/listener pairs and ng is the number of 
blocks, i. e. the number of times each talker/listener 
pair is used with each circuit. It is clear that 
the smallest possible 95% confidence interval is 
-|- 1.96 y V(C)/nc when nz is indefinitely large. There 
is, however, little advantage in increasing ny 
beyond the point where V,/ncnp= V(C)/nc because 
this reduces the confidence interval to within a 
factor of two of its minimum value. Greater preci- 
sion can more easily be achieved by increasing nc. 
The necessary number of observations can, and 
the necessary size of crew can, of course, only be 
decided after some measurements have been made 
and the results analysed. When a large experiment 
is to be undertaken much wasted effort can be avoid- 
ed by conducting a small-scale pilot experiment 


Table VITI 


j 
experiment | 


no. H 

10 We9)  |38.77 Ope? Sit 1.0 0 
1L* oiG- i. 0 Olwiel 87 ‘ae PTET 
12 (3.6 = Bee tS 0 
15 | 9.4 = = 0 as 0 
18 2.8 £2 = Os40u 0 
20 fl 600i = eae 20:90) = 0 
22 | 3.5 ee he 0 
my AY eee aS Ry eee 0 
24 3.0 ae . DAG = 0 


0 indicates that the component was not significant, 
— indicates that the design did not permit the com- 
ponent to be separated, 
* continuous spectrum noise was used instead of 
speech, 


to determine the order of magnitude of the com- 
ponents of variance and then arranging the main 
experiment accordingly. Similarly it is often not 
profitable to spend effort in improving the precision 
of each elementary value of m, for example by in- 
creasing the number of decisions per set. 

Table VIII gives a summary of the components 
of variance obtained in some recent experiments 
covering a wide range of telephone channels. Hight 
settings of sensitivity at 2 db intervals were used 
except for Nos. 22, 23 and 24 when there were six 
settings at 3 dh intervals. 


5. Conclusions 


The examples given here illustrate the use of meth- 
ods of conducting subjective acoustical measure- 
ments which are believed to be novel in this field. 
The methods are, however, developed from those 
in general use in biological experimentation. The 
methods depend on: 1. the use of appropriate ma- 
thematical models, 2. a systematic arrangement of 
the factors which are being controlled, and 3. ran- 
domization of the observations so that all uncon- 
trolled influences operate at random. Much waste 
of effort can be avoided by efficient design and ap- 
propriate analysis of the experiment. 
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BEITRAGE ZUR* ULT RASCHALLMESSTECHNIK IN FLUSSIGKEITEN 


Von JAN KOPPELMANN_ 
IIL. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 

Im ersten Teil der vorliegenden Arbeit wird ein neuartiges Sondenmikrophon zur punktfor- 
migen Abtastung von Ultraschallfeldern in Flissigkeiten bei Frequenzen oberhalb 100 kHz be- 
schrieben. 

Der zweite Teil befa®t sich im Anschlu8 an die Arbeit von HERTZ und MENDE [1] mit dem 
Zustandekommen des LANGEVINschen Schallstrahlungsdruckes. Es wird gezeigt, daB sich beim 
Einschalten des Ultraschallgenerators die Dichte der Flissigkeit im Schallstrahl erhéht, indem 
von aufBen Flissigkeit in den Schallstrahl nachstromt. 


Summary 

The first part of this paper describes a new form of microphone for point-to-point probing of 
ultrasonic fields in liquids at frequencies above 100 ke/s. 

The second part is concerned, as an extension to the previous work of HERTZ and MENDE [1], 
with the conditions for LANGEVIN’s formula for acoustic radiation pressure. It is shown that the 
introduction of the ultrasonic generator causes the density of the liquid in the sound beam to rise, 


following on the inflow of liquid from outside into the sound beam. 


Sommaire 


Dans la premiére partie du présent travail, on décrit un microphone-sonde d’un nouveau amo- 
déle pour l’étude, point par point, du champ ultra-sonore dans les liquides, aux fréquences su- 


périeures & 100 kHz. 


Dans la seconde partie, on s’occupe en complément du travail de HERTZ et MENDE [1], de 
l’établissement de la pression de radiation sonore de LANGEVIN. On montre que lorsque l’on met 
en route le générateur d’ultra-sons, la densité du liquide augmente dans le rayon sonore, du fait 
d’un transport de liquide de la périphérie a l’axe du pinceau sonore. 


1. Sondenmikrophone fiir Ultraschall 


Zur punktformigen Ausmessung von Ultraschall- 
feldern benétigt man ein Mikrophon, dessen Ab- 
messungen moglichst klein im Verhaltnis zur Wel- 
lenlinge sind. Abb. 1 zeigt ein Sondenmikrophon 
zur Messung von Ultraschallfeldern in Flissig- 
keiten bei Frequenzen oberhalb 100 kHz, das 
diesen Anforderungen gentigt. Ein etwa 0,5 mm 
dicker, 200 mm langer Stahldraht ist mit einem 
Gummi- oder Kunststoffschlauch tiberzogen, der 
den Schall méglichst schlecht leitet. Zwischen 
Draht und Schlauch befindet sich eine diinne, 
kérperschallisolierende Luftschicht. Am unteren 
Ende sind Draht und Schlauch zusammengekittet, 
um ein Eindringen von Flissigkeit zu verhindern. 

Wird die Sonde senkrecht in die Fliissigkeit ge- 
taucht, so wird die Spitze des Drahtes vom Schall- 
feld zu Schwingungen angeregt, die durch den 
Draht auf einen auBerhalb der Flissigkeit am 
oberen Ende des Drahtes angebrachten piezoelektri- 
schen Kristall tibertragen werden. Die erzeugte 
Wechselspannung wird in der tiblichen Weise ver- 
stirkt und mit einem Réhrenvoltmeter gemessen. 


Die Schallisolation zwischen dem Gummischlauch, 
der sich ebenfalls im Schallfeld befindet, und der 
Drahtspitze laBt sich in einfacher Weise tiberprii- 
fen,indem man die waagerecht liegende Sonde durch 
einen spitzen Ultraschallsprudel fahrt (siehe Abb. 2). 


schallempfindliche Spitze Schlauch 
Abb. 1. 


Schematischer Schnitt durch das Sondenmikrophon. 


Sondenmikrophon 

Abb. 2. 
Zur Uberpriifung der Schallisolation zwi- 
schen Schlauch und Drahtspitze wird das 
waagerecht liegende Mikrophon durch 
einen spitzen Ultraschallsprudel gefahren. 


Bei einem Sondenmikrophon mit einem Stahl- 
draht von 0,5 mm und einem Schlauch von 2... 3 
mm Durchmesser sinkt bei 500 kHz die vom Mikro- 
phon abgegebene Spannung auf den zehnten Teil, 
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wenn der Schlauch in mehr als 2... 3 mm Abstand 
von der Spitze vom Schall getroffen wird. Steht 
das Mikrophon senkrecht zum Schallstrahl, so zeigt 
es bei derselben Frequenz (3mm Wellenlange in 
Wasser) in einer stehenden Welle ein Amplituden- 
verhaltnis von rund 1: 20 an. Das Mikrophon ar- 
beitet in einem Bereich von drei Zehnerpotenzen 
des Schallwechseldruckes praktisch linear. Unter Be- 
nutzung eines guten Resonanzverstarkers lassen sich 
Schalleistungen von 2,5 x 10-2 W/cm? noch bequem 
nachweisen. Infolge der grofen Empfindlichkeit 
lassen sich auch auBerhalb der Resonanzfrequenz 
eines Ultraschallquarzes Messungen durchfihren. 


Abb, 3. Blick auf das mit Wasser gefiillte MeBgefaB. GroBe: 
30 x 35 x 45cm’. Vor den Glaswanden stehen 
Absorber aus schraggestellten Gummiplatten, um 
die Ausbildung stehender Wellen zu verhindern. 
Die Absorber schlucken etwa 97°/, der einfallenden 
Ultraschallenergie (bei 500 kHz).-In der Mitte der 
linken Wand sitzt der im Bild nicht sichtbare 
Ultraschallquarz. Mit dem Sondenmikrophon, 
dessen Spitze in das Wasser eintaucht, wird das 
Schallfeld ausgemessen. 


Abb. 3 zeigt eine Anordnung zur Ausmessung 
von Ultraschallfeldern. Mit dieser Anordnung wurde 
z. B. das Schallfeld eines Quarzes von 40 mm Durch- 
messer mit einer Resonanzfrequenz von 500 kHz 
auBerhalb der Resonanzfrequenz bei 200 kHz aus- 
gemessen. Das Ergebnis ist in Abb. 4 niedergelegt. 

Wahrend das Mikrophon fiir Relativmessungen 
in Schallfeldern bei fester Frequenz recht geeignet 
erscheint, ist im Hinblick auf die Absoluteichung 
zu beachten, daB sich auf dem Stahldraht stehende 
Wellen ausbilden und da8 sowohl Biegewellen wie 
auch Dehnwellen angeregt werden kénnen. Man darf 
daher keinen glatten Frequenzverlauf der Mikro- 
phonempfindlichkeit erwarten. AuBerdem 1a8t sich 
im Felde stehender Wellen leicht zeigen, daB je 


nach Frequenz, geometrischen Dimensionen der 


Stahlsonde und der Anordnung des Piezokristalles 
das Mikrophon in einem Schnellebauch, einem 
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a=304¢cm 


a=25,4cm 


a=204cm 


= 104cm 


Grosse des Quarzes | ern} 


Abb. 4. Schallwechseldruckverteilung senkrecht zumSchall- 
strahl bei verschiedenen Absténden vom Quarz. 
a= Abstand vom Quarz, r = senkrechter Abstand 
von der Mittelachse des Schallstrahles, p = Schall- 
wechseldruck in willktirlichen Einheiten, o = Me8- 
punkte. Auf den ersten beiden Kurven macht sich 
das Nahfeld des Quarzes bemerkbar (f = 200 kHz). 


Druckbauch oder in einem dazwischenliegenden 
Punkt seine maximale Empfindlichkeit erreicht. 

Dieser Ubelstand 148t sich bis zu einem gewissen 
Grade beheben, wenn man den Stahldraht durch 
einen Nickeldraht ersetzt. Infolge des Magneto- 
striktionseffektes wird in einer kleinen, um den 
Draht gelegten Induktionsspule eine der Ampli- 
tude der Dehnwelle auf dem Draht proportionale 
Wechselspannung erzeugt (siehe Abb. 5). Die Dehn- 
welle wird bevorzugt in den Druckbauchen ange- 
regt, so da8 solche Nickelsonden als saubere Druck- 
empfanger anzusprechen sind. Ferner lassen sich 


Nickeldraht 


Spule 


Abb. 5. Schematischer Schnitt durch ein Sondenmikrophon 
mit einem magnetostriktiven Nickeldraht. 
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die stehenden Wellen auf dem Draht weitgehend 
beseitigen, indem man das freie Ende des Drahtes 
mit Klebwachs oder dergleichen akustisch dampft. 

Der Restmagnetismus im Nickel reicht fiir die 
Benutzung des Nickeldrahtes im Sondenmikrophon 
aus. Bei exakten Messungen ist esratsam, die Magne- 
tisierung durch ein duBeres Magnetfeld konstant zu 
halten, da die Méglichkeit besteht, daB infolge Ultra- 
schalleinwirkung eine Entmagnetisierung eintritt. 


Abb. 6. Schematischer Schnitt durch einen Druckgra- 
dientenempfanger. Die beiden Nickeldraihte wer- 
den dicht nebenecinander in einem Schlauch unter- 
gebracht. 


Zwei Nickelsonden lassen sich zu einem Druck- 
gradientenempfanger zusammenbauen (s. Abb. 6). 
Die Induktionsspulen auf beiden Drahten sind 
hintereinandergeschaltet. Da sich auf den beiden 
Nickeldrahten stehende Wellen ausbilden, lassen 
sich durch Verschieben der Spulen Betrag und 
Phase der induzierten Spannung so abgleichen, 
daB die resultierende, vom Mikrophon abgegebene 
Spannung verschwindet. Nimmt man den Null- 
abgleich vor, wenn die Verbindungslinie der beiden 
Drahtspitzen senkrecht zum Schallstrahl steht, so 
erhalt man die fiir einen Gradientenempfanger ty- 
pische Richtcharakteristik. Abb. 7 zeigt ein Mef- 
ergebnis bei 500 kHz. 


Abb. 7. Richtcharakteristik des Druckgradientenempfin- 
gers bei 500 kHz in Wasser. Mittelpunktsabstand 
der Drahtspitzen 0,9 mm. In radialer Richtung ist 
die vom Mikrophon abgegebene Wechselspannung 
in willkilichen Einheiten aufgetragen. Die Un- 
symmetrie der Figur ist durch die Unsymmetrie 
der Dicke der Kittschicht zwischen Schlauch und 
Nickeldraéhten zu erklaren. 


Die Photographie Abb. 8 zeigt die bei den Mes- 
sungen benutzten Sondenmikrophone. Die Ent- 
wicklung der Mikrophone ist noch nicht abge- 


schlossen. Es sollte an dieser Stelle nur die Még- | 
lichkeit erortert werden, auf dem oben beschrie-— 
benen Wege Ultraschallfelder punktformig aus- 

zuumessen. ' 
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Abb. 8. Die bei den Messungen benutzten Mikrophone. 

Darunter liegt ein ZentimetermaBstab. 

1. Komplettes, piezoelektrisches Sondenmikrophon 
mit auswechselbarer Sonde, Drahtdurchmesser: 
1mm; 

. auswechselbare Sonde, Drahtdurchmesser: 
0,6 mm; 
3. auswechselbare Sonde, Drahtdurchmesser: 
0,3 mm; 
4. Magnetostriktive Sonde mit Spule, Drahtdurch- 
messer: 0,6 mm; 
. Magnetostriktive Sonde mit Spule, Drahtdurch- 
messer: 1mm; 
6. Druckgradientenmikrophon, Drahtdurchmesser 
0,6 mm, Mittelpunktsabstand der Drahte 0,9mm. 


bo 


or 


2. Zum Langevinschen Schallstrahlungsdruck 


Herrz und Menpe [1] haben gezeigt, daf der 
Lancervinsche Schallstrahlungsdruck nicht von der 
Richtung des Schallstrahles, sondern nur von 
dem Energiedichteunterschied der Energiedichte 
vor und hinter einer Grenzflache abhingt. Der 
Strahlungsdruck kommt dadurch zustande, dab 
sich die Dichte der Flissigkeit im Schallstrahl er- 
hoht, indem beim Eingchalten des Ultraschall- 
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generators Flissigkeit in den Schallstrahl stromt. 
Folgendes Experiment soll diese Auffassung be- 
statigen. 

Da der Schallstrahlungsdruck ein von GréoBen 
zweiter Ordnung abhingiger Effekt ist, 1iBt sich 
die dem Strahlungsdruck entsprechende Dichte- 
anderung statisch kaum nachweisen. Wird der 
Schallsender hingegen etwa 1000mal in der Se- 
kunde ein- und ausgeschaltet, so wird beim Kin- 
schalten Flissigkeit in den Schallstrahl hineinge- 
zogen. Beim Ausschalten erfolgt der Druckaus- 
gleich, indem dieselbe Wassermenge zurickflieBt. 
Das periodische Ein- und Ausschalten 1aft sich am 
einfachsten durch eine Amplitudenmodulation der 
am Ultraschallquarz liegenden Hochfrequenzspan- 
nung realisieren. In diesem Fall kann man die Be- 
grenzungsfliche des Schallstrahles als einen mit 
der Frequenz 1000 Hz radial schwingenden Zy- 
lindermantel auffassen. Der modulierte Schall- 
strahl strahlt also eine niederfrequente Schallwelle 
aus, die sich mit einem Mikrophon auBerhalb des 
hochfrequenten Schallstrahles bequem nachweisen 
14Bt. Der Schallwechseldruck ist allerdings um meh- 
rere GroBenordnungen kleiner als der Schallwech- 
seldruck des ,hochfrequenten Schallstrahles, da, 
wie bereits erwahnt, der Schallstrahlungsdruck ein 
Effekt zweiter Ordnung ist. Die Abstrahlung einer 
niederfrequenten Schallwelle ist ein eindeutiger 
Beweis fiir das Nachstromen der Fliissigkeit in den 
Schallstrahl. 

FaBt man den modulierten Schallstrahl als 
schwingenden Zylindermantel auf, so ist zunachst 
zu beachten, daB sich der Schall mit endlicher Ge- 
schwindigkeit ausbreitet. Demzufolge schwingt der 
Zylindermantel nicht auf der gesamten Lange 
gleichphasig, sondern die Phasenlage dndert sich 
mit der zu 1000 Hz gehorigen Wellenlange (in Was- 
ser 150 cm), siehe Abb. 9. Hat der Schallstrahl eine 
hinreichende Lange, so bedeutet das, daB die nieder- 
frequente Welle gerichtet, in Richtung des Schall- 
strahles, abgestrahlt wird. Infolge der groBeren 
Wellenlange dieser Schallwelle reicht das Nahfeld 
weit iiber den Ultraschallstrah! hinaus. 

Ferner fiihrt das Bild eines schwingenden Zylin- 
dermantels zu der Konsequenz, daB zu einem festen 
Zeitpunkt in einem festen Abstand vom Ultra- 
schallsender ein Uberdruck im Schallstrahl einem 
Unterdruck auBerhalb des Schallstrahles entspre- 
chen mu8 und umgekehrt. Halt man ein kleines, 
nur auf die Niederfrequenz ansprechendes Mikro- 
phon mit einer senkrecht zur Schallfortpflanzungs- 


_ richtung stehenden, druckempfindlichen, ebenen 


Flache in den Schallstrahl, so tritt beim Einschalten 
des Schallsenders an dieser Flache ein Uberdruck, 
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Abb. 9. Der modulierte Ultraschallstrahl laBt sich beztig- 
lich der Modulationsfrequenz als schwingender 
Zylindermantel auffassen. 


der Schallstrahlungsdruck, auf. Dieser kommt da- 
durch zustande, da Fliissigkeit in den Schallstrahl 
hineingezogen wird [1]. Befindet sich das Mikro- 
phon auBerhalb des Schallstrahles, so tritt dort 
zur gleichen Zeit notgedrungen ein Unterdruck auf. 
Wird das Mikrophon in einer Richtung senkrecht 
zum Schallstrahl bewegt, so beobachtet man an der 
Grenze des Ultraschallstrahles einen Phasensprung 
von 180°, der sich meBtechnisch leicht nachweisen 
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Abb. 10. Niederfrequenter Wechseldruck als Funktion des 
Abstandes r von der Mittelachse des Schallstrah- 
les im Fernfeld eines 500 kHz-Quarzes von 4cm 
Durchmesser. Die Modulationsfrequenz betragt 
etwa 1500 Hz. Der Phasensprung von 180° wird 
durch das Auftragen negativer Werte angedeutet ; 
o = MeBpunkte. 


148t. Abb. 10 zeigt das Ergebnis einer solchen Mes- 
sung. Bei tieferen Modulationsfrequenzen sinkt die 
Intensitat des niederfrequenten Wechseldruckes 
auBerhalb des Ultraschallstrahles stark ab, weil der 
Durchmesser des Schallstrahles und die GefaBdi- 
mensionen im Verhiltnis zur Wellenlainge zu klein 
sind. Bei héheren Frequenzen ist das gesamte MeBge- 
fiB mit niederfrequenten, stehenden Wellen erfillt. 
Die Anregung zu diesen Arbeiten ging von Herrn 
Prof. E. Mrver aus, dem ich an dieser Stelle fir 
viele, wertvolle Hinweise meinen besonderen Dank 

aussprechen mdéchte. 
(Eingegangen am 1. Juli 1951.) 

Schrifttum 

[1] HERTZ, G. und MENDE, H., Z. Phys. 114 [1939], 354. 
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International Congress on Electro-acoustics 
and 
International Symposium on the Sound Insulation 
of Light-weight Structures 
15—21 June 1953, Netherlands 


The Congress and Symposium are organized by the 
NETHERLANDS ACOUSTICAL SOCIETY (Geluidstichting) 
under the auspices of the International Commission on 
Acoustics (ICA, see Acustica 1 [1951], 143) of the Inter- 
national Union of Pure and Applied Physics. 


Organizing Committee: 
R. VERMEULEN, Eindhoven, chairman 
C. W. KOSTEN, Delft, vice-chairman 
W. Th. BAHLER, Delft 
J. J. GELUK, Hilversum 
H. MOL, Leiden 
P. A. de LANGE, Delft, secretary 


Congress-centre: Delft (near the Hague). 


Meetings: in Delft (Technical University), Hindhoven 
(Philips Works) and Hilversum (Broadcasting centre). 


Opening lecture by the President of the International 
Commission on Acoustics, R. H. BOLT, Cambridge, 
Mass. 


Sections: Apart from the general sessions members will 
meet in specialized sections. In each section a general 
review will be given and papers on original work, each 
lasting 15 minutes, read. The symposium on sound in- 
sulation is for organizational reasons listed as section 
VI of the Congress. A section on electro-acoustics 
applied to musical instruments will be added, if 
sufficient papers are forthcoming; general reviewer 
E. G. RICHARDSON, Newcastle-upon-Tyne. 


Nr. Subject matter of section General reviewer 


R. VERMEULEN, 

Eindhoven 
E. MEYER, G6ttingen 
L. L. BERANEK, Cam- 
bridge, Mass. 

P. CHAVASSE, Paris 
G. BRADFIELD, Ted- 
dington 


I Phonography and recording 
(dise, film, wire, tape) 

II Public address systems 

III Acoustic measurements 


IV Hearing aids and audiometers 
V_ Electroacoustics in ultrasonics 


VI Sound insulation of light- 


weight structures C. W. KOSTEN, Delft 


Papers: 
Those intending to read a paper are requested to com- 
municate as soon as possible name, address, section 
and title of their paper to the Congress secretary Mr 
P. A. de LANGE, Laboratorium voor Technische Phy- 
sica, Mijnbouwplein 11, Delft, Netherlands. The Com- 
mittee may have to limit the number of papers. 
Languages: 
English, French and German are the official languages. 
Membership: 
The fee for ordinary membership will be about 20—40 
guilders. The exact amount will be fixed later; it de- 
pends upon whether the papers will be printed or not. 
Fee for associate members: 15 guilders. 
Excursions, ete.: 
A number of technical and non-technical excursions 
are being prepared and special arrangements made for 
ladies accompanying members. 


Details: 


A circular giving details will be distributed in due time. 
Further announcements will be made in future issues 
of this Journal. 


For the Organizing Committee 


C. W. Kosten 
P. A. de Lange 


Cinquiéme Congrés International d’Oto- Rhino- La- 
ryngo- Broncho- Oesophagologie 
8—15 Juin, 1953, Amsterdam 


Organisé par la Société Néerlandaise d’Oto-Rhino-La- 
ryngologie. 

Prof. EELCO HUIZINGA, président 

W. H. STRUBEN, secrétaire général 

J. BIJTEL et P. G. GERLINGS, secrétaires scientifiques 

Miles A.C. SCHIPPERS et A. DOKKUM, secrétaires de 

réception 

Mme F. VELLEMAN-PINTO, trésoriére 
A coté des Réunions Combinées (sujets choisis) les Membres 
du Congrés seront invités 4 des Réunions Restreintes. Les 
sujets choisis sont: 


TUMEURS DES BRONCHES 


A. SOULAS et P. MOUNIER-KUHN, Tumeurs bénignes 
des bronches 

L. PIETRANTONI, La clinique et l’endoscopie des cancers 
des bronches 

L. H. CLERF, Cytology of tumours of the bronchus 


AUDIOLOGIE 


C. 8. HALLPIKE, Physiologie problems of Audiology 
G. E. SHAMBAUGH, Clinics of Audiology 
S. R. SILVERMAN, Social aspects of Audiology 


ALLERGIE 


G. DOHLMAN, Theory of Allergy 
F. H. HANSEL, Methods of Allergic Examination 
R. MELCHIOR, L’hypersensibilité d’origine nonallergique 


Programme des Dames. Un Comité de Dames s’occu- 
pera des dames étrangéres pendant leur séjour. 


Langues officielles: l’anglais, le frangais, lallemand. 


Communications: deux sortes de communications sont 
demandeées. 

a. Communications originales (15 minutes). Le nom- 
bre des communications originales est limité, ainsi que 
leur discussion. On est prié de s’adresser au Secrétaire 
(voir plus bas) pour les informations 4 ce sujet. 

b. Contributions aux discussions des sujets choi- 
sis ‘des Réunions Combinées: 5 minutes par 
discours, une heure par sujet choisi. 


Cotisations: 100 florins pour les membres, 25 florins pour 
les Dames et les invités. 


Logement & Amsterdam ou aux environs. 


Renseignements. Pour inscriptions, informations, sug- 
gestions, etc. on est prié de s’adresser au Secrétaire du 
Congrés. 

W. H. Struben, 
J.J. Viottastraat 1, 
Amsterdam-Z. 

(abrégé par la Rédaction) 


